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1. Introduction 

1.1 Overview. 

In the recent papers [1, 2], we presented a general theory of partition function zeros in models 
with periodic boundary conditions and interaction depending on one complex parameter. The 
analysis was based on a set of assumptions, called Assumptions A and B in [2], which are es- 
sentially statements concerning differentiability properties of certain free energies supplemented 
by appropriate non-degeneracy conditions. On the basis of these assumptions we characterized 
the topology of the resulting phase diagram and showed that the partition function zeros are in 
one-to-one correspondence with the solutions to specific (and simple) equations. In addition, the 
maximal degeneracy of the zeros was proved to be bounded by the number of thermodynamically 
stable phases, and the distance between the zeros and the corresponding solutions was shown to 
be generically exponentially small in the linear size of the system. 

The reliance on Assumptions A and B in [2] permitted us to split the analysis of partition 
function zeros into two parts, which are distinct in both mathematical and physical content: one 
concerning the zeros of a complex (in fact, analytic) function — namely the partition function 
with periodic boundary conditions — subject to specific requirements, and the other concerning 
the control of the partition function in a statistical mechanical model depending on one complex 
parameter. The former part of the analysis was carried out in [2]; the latter is the subject of this 
paper. Explicitly, the principal goal of this paper can be summarized as follows: We will define a 
large class of lattice spin models (which includes several well-known systems, e.g., the Ising and 
Blume-Capel models) and show that Assumptions A and B are satisfied for every model in this 
class. On the basis of [2], for any model in this class we then have complete control of the zeros 
of the partition function with periodic boundary conditions. 

The models we consider are characterized by two properties: the existence of only a finite 
number of ground states and the availability of a contour representation. In our setting, the term 
ground state will simply mean a constant — or, after some reinterpretations, a periodic — infinite 
volume spin configuration. Roughly speaking, the contour representation will be such that the 
contours correspond to finite, connected subsets of the lattice where the spin configuration differs 
from any of the possible ground states. A precise definition of these notions is a bit technical; 
details will be provided in Section 3. Besides these properties, there will also be a few quantitative 
requirements on the ground state energies and the scaling of the excess contour energy with the 
size of the contour — the Peierls condition — see Sections 2.1 and 3.2. 

These two characteristic properties enable us to apply Pirogov-Sinai theory — a general method 
for determining low-temperature properties of a statistical mechanical model by perturbing about 
zero-temperature. The first formulation of this perturbation technique [16, 17] applied to a class 
of models with real, positive weights. The original "Banach space" approach of [16, 17] was 
later replaced by inductive methods [9], which resulted in a complete classification of translation- 
invariant Gibbs states [21]. The inductive techniques also permitted a generalization of the char- 
acterization of phase stability/coexistence to models with complex weights [5]. However, most 
relevant for our purposes are the results of [6], dealing with finite-size scaling in the vicinity of 
first-order phase transitions. There Pirogov-Sinai theory was used to derive detailed asymptotics 
of finite volume partition functions. The present paper provides, among other things, a variant of 
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[6] that ensures appropriate differentiability of the so-called metastable free energies as required 
for the analysis of partition function zeros. 

The remainder of this paper is organized as follows. Section 1.2 outlines the class of models 
of interest. Section 1.3 defines the ground state and excitation energies and introduces the torus 
partition function — the main object of interest in this paper. Section 2.1 lists the assumptions 
on the models and Section 2.2 gives the statements of the main results of this paper. These 
immediately imply Assumptions A and B of [2] for all models in the class considered. Sections 3 
and 4 introduce the necessary tools from Pirogov-Sinai theory. These are applied in Section 5 to 
prove the main results of the paper. 

1.2 Models of interest. 

Here we define the class of models to be considered in this paper. Most of what is to follow in 
this and the forthcoming sections is inspired by classic texts on spin models, Gibbs states and 
Pirogov-Sinai theory, e.g., [8, 18, 20, 21]. 

We will consider finite-state spin models on the d-dimensional hypercubic lattice Z d for d > 2. 
At each site x£Z' i the spin, denoted by a x , will take values in a finite set S. A spin configuration 
a = (o~ x ) x( z Z d is an assignment of a spin to each site of the lattice. The interaction Hamiltonian 
will be described using a collection of potentials (3>a), where A runs over all finite subsets of Z d . 
The $a are functions on configurations from 5 Z with the following properties: 

(1) The value $a(o~) depends only on a x with x G A. 

(2) The potential is translation invariant, i.e., if o' is a translate of a and A' is the correspond- 
ing translate of A, then &a'(o~) = ^a(c')- 

(3) There exists an R > 1 such that <I>a = for all A with diameter exceeding R+l. 

Here the diameter of a cubic box with L x • • • x L sites is defined to be L while for a general A C 
7L d it is the diameter of the smallest cubic box containing A. The constant R is called the range 
of the interaction. 

Remark 1. Condition (2) has been included mostly for convenience of exposition. In fact, all of 
the results of this paper hold under the assumption that $a are periodic in the sense that <I>a' (<t) = 
<&a(c') holds for A and a related to A' and a' by a translation from {aL) d for some fixed integer a. 
This is seen by noting that the periodic cases can always be converted to translation-invariant ones 
by considering block-spin variables and integrated potentials. 

As usual, the energy of a spin configuration is specified by the Hamiltonian. Formally, the 
Hamiltonian is represented by a collection of functions (f3H\) indexed by finite subsets of Z d , 
where (3H\ is defined by the formula 



(The superfluous f3, playing the role of the inverse temperature, appears only to maintain formal 
correspondence with the fundamental formulas of statistical mechanics.) In light of our restriction 
to finite-range interactions, the sum is always finite. 

We proceed by listing a few well known examples of models in the above class. With the 
exception of the second example, the range of each interaction is equal to 1: 





(1.1) 



A' : A'nA^0 
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Ising model. Here S = { — 1,+1} and <J>a(c) ^ only for A containing a single site or a 
nearest-neighbor pair. In this case we have 

*a{*) = l'*"" if t = | X \ , , (1-2) 

y-JaxO-y, if A = {x,y\ with \x - y\ = 1. 

Here J is the coupling constant, h is an external field and \x — y\ denotes the Euclidean distance 
between x and y. 

Perturbed Ising model. Again S = {— 1, +1}, but now we allow for arbitrary finite range 
perturbations. Explicitly, 

<& A (a) = {-^' ifA = W, 

\-JaY\x G k°x if |A| > 2and diamA< R + l. 

The coupling constants J\ are assumed to be translation invariant (i.e., J\ = J\/ if A and A' are 
translates of each other). The constant h is again the external field. 

Blume-Capel model. In this case S = { — 1, 0, +1} and &a(o~) = unless A is just a single site 
or a nearest-neighbor pair. Explicitly, we have 

-Xa^ — ha x , ifA = {x}, 

J{o- x ~ o-y) 2 , if A = {x,y} with \x - y\ = 1. 

Here J is the coupling constant, A is a parameter favoring ±1 against 0-spins and h is an external 
field splitting the symmetry between +1 and —1. 

Potts model in an external field. The state space has q elements, S = {1, . . . , q} and <1>a is again 
nontrivial only if A is a one-element set or a pair of nearest-neighbor sites. Explicitly, 

*a(<t) = <[" W '" 1 ' ifA = W' (1.5) 

{-JSa x ,a y , if A = {x,y} with|x-y| = 1. 

Here S a y equals one if a = a' and zero otherwise, J is the coupling constant and h is an 
external field favoring spin value 1. Actually, the results of this paper will hold only for the low- 
temperature regime (which in our parametrization corresponds to J » loggO; a more general 
argument covering all temperatures (but under the condition that q is sufficiently large) will be 
presented elsewhere [3, 4]. 

Any of the constants appearing in the above Hamiltonian can in principle be complex. How- 
ever, not all complex values of, e.g., the coupling constant will be permitted by our additional 
restrictions. See Section 2.3 for more discussion. 

1.3 Ground states, excitations and torus partition function. 

The key idea underlying our formulation is that constant configurations represent the potential 
ground states of the system. (A precise statement of this fact appears in Assumption C2 below.) 
This motivates us to define the dimensionless ground state energy density e m associated with spin 
m G S by the formula 

e ™= E m* A(cjm) ' (L6) 

A: A30 ' ' 
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where | A| denotes the cardinality of the set A and where a m is the spin configuration that is equal 
to m at every site. By our restriction to finite-range interactions, the sum is effectively finite. 

The constant configurations represent the states with minimal energy; all other configurations 
are to be regarded as excitations. Given a spin configuration a, let -Br(ct) denote the union of all 
cubic boxes A C 2, d of diameter 2R + 1 such that a is not constant in A. We think of Br((j) as 
the set on which a is "bad" in the sense that it is not a ground state at scale R. The set -Br(cx) 
will be referred to as the R-boundary of a. Then the excitation energy E(a) of configuration a 
is defined by 

E (*) = E E jxi*A(*). a- 7 ) 

xeB R (a)A: xeA 

To ensure that the sum is finite (and therefore meaningful) we will only consider the configura- 
tions a for which Br(ct) is a finite set. 

The main quantity of interest in this paper is the partition function with periodic boundary 
conditions which we now define. Let L > 2R + 1, and let denote the torus of L x L x 
• • • x L sites in 'L d , which can be thought of as the factor of Z d with respect to the action of the 
subgroup (LZ) d . Let us consider the Hamiltonian (3Hl : S L — > C defined by 

(3H L (a)= $ a(<t), aeS TL , (1.8) 

A: AcT L 

where $a w& retractions of the corresponding potentials from Z rf to T^. (Here we use the trans- 
lation invariance of $a.) Then the partition function with periodic boundary conditions in is 
defined by 

Z p L er = ^ e~ pH ^ a \ (1.9) 

In general, Z^ er is a complex quantity which depends on all parameters of the Hamiltonian. 
We note that various other partition functions will play an important role throughout this paper. 
However, none of these will be needed for the statement of our main results in Section 2, so we 
postpone the additional definitions and discussion to Section 4. 

We conclude this section with a remark concerning the interchangeability of the various spin 
states. There are natural examples (e.g., the Potts model) where several spin values are virtually 
indistinguishable from each other. To express this property mathematically, we will consider the 
situation where there exists a subgroup of the permutations of S such that if it € <25 then 
e 7r(m) = e m and E(ir(a)) = E(a) for each m£5 and each configuration a with finite Br{g), 
where ir(cr) is the spin configuration taking value Tr(a x ) at each x. (Note that Br{k{o-)) = Br(ct) 
for any such permutation tt.) Then we call two spin states m and n interchangeable if m and n 
belong to the same orbit of the group (5 on S. 

While this extra symmetry has absolutely no effect on the contour analysis of the torus partition 
sum, it turns out that interchangeable spin states cannot be treated separately in our analysis of 
partition function zeros. (The precise reason is that interchangeable spin states would violate 
our non-degeneracy conditions; see Assumption C3-C4 and Theorem A3-4 below.) To avoid this 
difficulty, we will use the factor set 1Z = 5/0 instead of the original index set S when stating 
our assumptions and results. In accordance with the notation of [2], we will also use r to denote 
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the cardinality of the set 1Z, i.e., 1Z = {1, 2, . . . , r}, and q m to denote the cardinality of the orbit 
corresponding torn S 71. 

2. Assumptions and results 

In this section we list our precise assumptions on the models of interest and state the main results 
of this paper. 

2.1 Assumptions. 

We will consider the setup outlined in Sections 1.2-1.3 with the additional assumption that the 
parameters of the Hamiltonian depend on one complex parameter z which varies in some open 
subset G of the complex plane. Typically, we will take z = e h or z = e 2h where h is an external 
field; see the examples at the end of Section 1.2. Throughout this paper we will assume that the 
spin space S, the factor set 1Z, the integers q m and the range of the interaction are independent of 
the parameter z. We will also assume that the spatial dimension d is no less than two. 

The assumptions below will be expressed in terms of complex derivatives with respect to z. 
For brevity of exposition, let us use the standard notation 

0* = 3 and fc = 3(&+*&) t 2 - 1 ) 

for the derivatives with respect to z and z, respectively. Here x = ^tez and y = Qmz. Our 
assumptions will be formulated for the exponential weights 

<p A (a,z) = e-^ z \ p z (a)=e~ E ^ and 9 m (z) = e~^ z \ (2.2) 

where we have now made the dependence on z notationally explicit. In terms of the 6 m 's and the 
quantity 

9{z) = max|Mz)| (2-3) 

we define the set 5£ a (m) by 

JS? Q (m) = {z£0: \9 m (z)\ > 6{z)e a ). (2.4) 

Informally, J£ a (m) is the set of z for which m is "almost" a ground state of the Hamiltonian. 

Since we want to refer back to Assumptions A and B of [2], we will call our new hypothesis 
Assumption C. 

Assumption C. There exist a domain & C C and constants a,M,r G (0, oo) such that the 
following conditions are satisfied. 

(0) For each a G and each finite A C 1, d , the function z i— ► v?a(<7, z) is holomorphic in G. 

(1) For all m G S, all z G & and all t = 0, 1, 2, the ground state weights obey the bounds 

\d e z e m (z)\ < M e 6(z) (2.5) 

In addition, the quantity 0(z) is uniformly bounded away from zero in &. 

(2) For every configuration a with finite i?-boundary Br((j), the Peierls condition 

\d e zPz (a)\ < {M\B R (a)\Y{e-^(z)f R{ ^ (2.6) 
holds for all z G & and I = 0, 1, 2. 
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(3) For all distinct m, n G 1Z and all z G =Sf a (m) n Jf a (n), we have 

d z e m {z) d z e n (z 



> a. (2.7) 



inf< 



9 m (z) 9 n (z) 

(4) If Q C 7?. is such that \Q\ > 3, then for any z G f) m( zQ S£ a (m) we assume that the complex 
quantities v m {z) = 6 m {z)~ l d z 9 m {z), m G Q, regarded as vectors in M 2 , are vertices of a 
strictly convex polygon. Explicitly, we demand that the bound 

■ < v m (z) - ^2 u n v n {z) : uj h > 0, ^ w„ = 1 > > a (2.8) 

neQ\{m} rteQ\{m} 

holds for every meQ and every z G PlneC -^a( n )- 

Assumptions CO-2 are very natural; indeed, they are typically a consequence of the fact that the 
potentials <^a(<7, z) — and hence also 9 m (z) and p z (a) — arise by analytic continuation from the 
positive real axis. Assumptions C3-4 replace the "standard" multidimensional non-degeneracy 
conditions which are typically introduced to control the topological structure of the phase dia- 
gram, see e.g. [16, 17, 20]. (However, unlike for the "standard" non-degeneracy conditions, here 
this control requires a good deal of extra work, see [2].) Assumption C4 is only important in the 
vicinity of multiple coexistence points (see Section 3.2); otherwise, it can be omitted. 

Remark 2. For many models, including the first three of our examples, the partition function has 
both zeros and poles, and sometimes even involves non-integer powers of z. In this situation it is 
convenient to multiply the partition function by a suitable power of z to obtain a function that is 
analytic in a larger domain. Typically, this different normalization also leads to a larger domain & 
for which Assumption C holds. Taking, e.g., the Ising model with z = e 2h , one easily verifies that 
for low enough temperatures, Assumption C holds everywhere in the complex plane — provided 
we replace the term — ha x by —h(a x + 1). By contrast, in the original representation (where 
<P{ x }{cr, z) = (yfz~y x ), one needs to take out a neighborhood of the negative real axis (or any 
other ray from zero to infinity) to achieve the analyticity required by Assumption CO. 

Remark 3. If we replace the term — ha x in (1.2-1.4) by — h(a x + l), Assumption C (with z = e 2h 
for the Ising models, and z = e h for the Blume Capel and Potts model) holds for all four examples 
listed in Section 1.2, provided that the nearest-neighbor couplings are ferromagnetic and the 
temperature is low enough. (For the perturbed Ising model, one also needs that the nearest- 
neighbor coupling is sufficiently dominant.) 



2.2 Main results. 

Now we are in a position to state our main results, which show that Assumptions A and B from [2] 
are satisfied and hence our conclusions concerning the partition function zeros hold. The structure 
of these theorems parallels the structure of Assumptions A and B. We caution the reader that the 
precise statement of these results is quite technical. For a discussion of the implications of these 
theorems, see Section 2.3. The first theorem establishes the existence of metastable free energies 
and their relation to the quantities 9 m . 
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Theorem A Let M G (0, oo) and a G (0, oo). Then there is a constant tq depending on M, a, 
the number of spin states \S\ and the dimension d such that if Assumption C holds for the con- 
stants M, a, some open domain ^CC and some r > tq, then there are functions Q m : G — > C, 
m G TZ,for which the following holds: 

(1) There are functions s m : & — > C, m G TZ, such that Q m {z) can be expressed as 

(m(z) = e m (z)e s ™W and \s m (z)\ < e"^ 2 . (2.9) 

In particular, the quantity Q(z) = max mG 7j | ( m (z)\ is uniformly positive in iff. 

(2) Each function ( m , viewed as a function of two real variables x = 3ftez and y = 9mz, is twice 
continuously differ entiable on and satisfies the Cauchy-Riemann equations dz( m (z) = 
for all z G ,5^ m , where 

y m = {ze0:\(m(z)\ = ((z)}. (2.10) 

In particular, ( m is analytic in the interior of 5^ m . 

(3) For any pair of distinct indices m, n G TZ and any z G 5^ m Pi ,5^ n we have 

dzCmjz) _ d z ( n {z) 
Crn(z) Cn(z) 

(4) If Q C TZ is such that \Q\ > 3, then for any z G f] m£ Q 

v m (z) = , m G Q, (2.12) 

(,m\Z) 

are the vertices of a strictly convex polygon in C ~ M 2 . 

Theorem A ensures the validity of Assumption A in [2] for any model satisfying Assumption C 
with r sufficiently large. Assumption A, in turn, allows us to establish several properties of the 
topology of the phase diagram, see Section 2.3 below for more details. 

Following [2], we will refer to the indices in TZ as phases, and call a phase m G TZ stable 
at z if |Cm(z)| = C{z). We will say that a point z G & is a point of phase coexistence if there 
are at least two phases m G TZ which are stable at z. In [2] we introduced these definitions 
without further motivation, anticipating, however, the present work which provides the technical 
justification of these concepts. Indeed, using the expansion techniques developed in Sections 3 
and 4, one can show that, for each m G S that corresponds to a stable phase in TZ, the finite 
volume states with m-boundary conditions tend to a unique infinite- volume limit (-} m in the 
sense of weak convergence on linear functional on local observables. (Here a local observable 
refers to a function depending only on a finite number of spins). The limit state is invariant 
under translations of Z d , exhibits exponential clustering, and is a small perturbation of the ground 
state a m in the sense that (5 ax> k)m = 5 m ,k + 0(er T l 2 ) for all x G Z d . 

Remark 4. Note that two states (-) m and (-) m i are considered as two different versions of the 
same phase if m and m! are indistinguishable, in accordance with our convention that TZ, and 
not S, labels phases. Accordingly, the term phase coexistence refers to the coexistence of dis- 
tinguishable phases, and not to the coexistence of two states labelled by different indices in the 
same orbit TZ. This interpretation of a "thermodynamic phase" agrees with that used in physics, 
but disagrees with that sometimes used in the mathematical physics literature. 



> a - 2e" 



-T/2 



(2.11) 
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While Theorem A is valid in the whole domain G, our next theorem will require that we restrict 
ourselves to a subset G C G with the property that there exists some e > such that for each 
point z £ G, the disc D e (z) of radius e centered at z is contained in G. (Note that this condition 
requires G to be a strict subset of G, unless G consists of the whole complex plane). In order 
to state the next theorem, we will need to recall some notation from [2]. Given any m £ 1Z and 
5 > 0, let ^(m) denote the region where the phase m is "almost stable," 

y s (m) = {z£G: \Cm(z)\ > e' 5 ((z)}. (2.13) 

For any Q C 1Z, we also introduce the region where all phases from Q are "almost stable" while 
the remaining ones are not, 

%{Q) = p| y 5 {m) \ |J y^nj, (2.14) 

meQ neQ c 

with the bar denoting the set closure. 

Theorem B Let M, a, e G (0, oo), and let r > tq, where tq is the constant from Theorem A, 
and let k = r/4. Let G C C and G <Z G be open domains such that that Assumption C holds 
in G and H> e (z) C G for all z £ G. Then there are constants Co (depending only on M), Mq 
(depending on M and e), and Lq (depending on d, M, r and e) such that for each m G 1Z 
and each L > Lq there is a function (f^ : ^/^(m) — > C such that the following holds for all 
L > L : 

(1) The function Z^ CT is analytic in G. 

(2) Each ci L) is non-vanishing and analytic in ^/^(m). Furthermore, 



and 



d z \og 



log 



Cm(z) 



< e 



-rL/8 



Cm (2:) 



+ 



<9 2 log 



ti\z) 



Cm{z) 



< e 



-tL/8 



(2.15) 



(2.16) 



hold for allm £ 1Z and all z £ ^ t /L( m -)- 
(3) For each m £ 1Z, all t > 1, and all z £ ^/^(m), we have 



dlti\z) 



tt\z) 



Moreover, for all distinct m,n £ 1Z and all z £ y K /L( m ) H ^ K /z,(n) 

>a-2e' r ' 2 . 



d z d L \z) 



&\z) 



d L \z) 



(4) For any Q C TZ, the difference 



meQ 



■,[tiXz)\ 



(2.17) 



(2.18) 



(2.19) 
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satisfies the bound 



Hs Q , L (z)\<iKCoL d Y +1 az) Ld ( J> m e 



tL/W 



(2.20) 



men 



for all£>0 and all z € W K / L {Q)- 

Theorem B proves the validity of Assumption B from [2]. Together with Theorem A, this in 
turn allows us to give a detailed description of the positions of the partition function zeros for all 
models in our class, see Section 2.3. 

The principal result of Theorem B is stated in part (4): The torus partition function can be 
approximated by a finite sum of terms — one for each "almost stable" phase m G 1Z — which have 
well controlled analyticity properties. As a consequence, the zeros of the partition function arise 
as a result of destructive interference between almost stable phases, and all zeros are near to the 
set of coexistence points = Um^n ^ see Section 2.3 for further details. Representa- 
tions of the form (2.19) were crucial for the analysis of finite-size scaling near first-order phase 
transitions [6]. The original derivation goes back to [5]. In our case the situation is complicated 
by the requirement of analyticity; hence the restriction to z G ^ k /l(Q) in (4)- 

2.3 Discussion. 

As mentioned previously, Theorems A and B imply the validity of Assumptions A and B of [2], 
which in turn imply the principal conclusions of [2] for any model of the kind introduced in 
Section 1.2 that satisfies Assumption C with r sufficiently large. Instead of giving the full state- 
ments of the results of [2], we will only describe these theorems on a qualitative level. Readers 
interested in more details are referred to Section 2 of [2]. 

Our first result concerns the set of coexistence points, & = {J m ^ n y m H .5^1, giving rise to the 
complex phase diagram. Here Theorem 2. 1 of [2] asserts that is the union of a set of simple, 
smooth (open and closed) curves such that exactly two phases coexist at any interior point of the 
curve, while at least three phases coexist at the endpoints — these are the multiple points. More- 
over, in each compact set, any two such curves cannot get too close without intersecting and there 
are only a finite number of multiple points. These properties are of course direct consequences of 
the non-degeneracy conditions expressed in Theorem A3-4. 

Having discussed the phase diagram, we can now turn our attention to the zeros of Z£ cr '. The 
combined results of Theorems 2.2-2.4 of [2] yield the following: First, all zeros lie within 0(L~ d ) 
of the set Second, along the two-phase coexistence lines with stable phases m, n G 1Z, the 
zeros are within 0(e~ cL ), for some c > 0, of the solutions to the equations 



Consecutive solutions to these equations are separated by distances of order L , i.e., there are of 
the other L d zeros per unit length of the coexistence line. Scaling by L d , this allows us to define 
a density of zeros along each two-phase coexistence line, which in the limit L —> oo turns out to 
be a smooth function varying only over distances of order one. 



qU Ld \Uz)\=q 1 J Ld \Uz)\, 
L d Arg(( m (z)/( n {z)) =7rmod27r. 



(2.21) 



(2.22) 
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Figure 1. A schematic figure of the solutions to (2.21-2.22) giving the approximate locations of 
partition function zeros of the Ising model in parameter z which is related to the external field h 
by z — e 2h . The plot corresponds to dimension d — 2 and torus side L = 8. The expansion used for 
calculating the quantities f± is shown in (2.23). To make the non-uniformity of the spacing between 
zeros more apparent, the plot has been rendered for the choice e 2J = 2.5 even though this is beyond 
the region where we can prove convergence of our expansions. 



Near the multiple points the zeros are still in one-to-one correspondence with the solutions of 
a certain equation. However, our control of the errors here is less precise than in the two-phase 
coexistence region. In any case, all zeros are at most (r — 1) -times degenerate. In addition, for 
models with an Ising-like plus-minus symmetry, Theorem 2.5 of [2] gives conditions under which 
zeros will lie exactly on the unit circle. This is the local Lee- Yang theorem. 

Let us demonstrate these results in the context of some of our examples from Section 1.2. We 
will begin with the standard Ising model at low temperatures. In this case there are two possible 
phases, labeled + and — , with the corresponding metas table free energies given as functions 

of z = e 2h by 

C±(z) = exp{±/i + e - 2dJ ^ 2h + 0(e-( 4d - 2 ) J )}. (2.23) 

Symmetry considerations now imply that |C+(^)| = IC-WI if an d only if 3te/i = 0, i.e., \z\ = 
1, and, as already known from the celebrated Lee- Yang Circle Theorem [11], the same is true 
for the actual zeros of Z^ 01 . However, our analysis allows us to go further and approximately 
calculate the solutions to the system (2.21-2.22), which shows that the zeros of Z^ er lie near the 
points z = e t9k , where k = 0, 1, . . . , L d — 1 and 



Ok — 



2k + 1 



7T + 2e~ 2dJ sin 



2k + 1 



-7T 



+ 0(e 



-(4<2-2)J 



)• 



(2.24) 



L d V L d 

Of course, as L increases, higher and higher-order terms in e~ J are needed to pinpoint the location 
of any particular zero (given that the distance of close zeros is of the order L~ d ). Thus, rather than 
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providing the precise location of any given zero, the above formula should be used to calculate 
the quantity 6u+i — Ok, which is essentially the distance between two consecutive zeros. The 
resulting derivation of the density of zeros is new even in the case of the standard Ising model. A 
qualitative picture of how the zeros span the unit circle is provided in Fig. 1 . 

A similar discussion applies to the "perturbed" Ising model, provided the nearest-neighbor 
coupling is ferromagnetic and the remaining terms in the Hamiltonian are small in some appro- 
priate norm. In the case of general multi-body couplings, the zeros will lie on a closed curve 
which, generically, is not a circle. (For instance, this is easily verified for the three-body inter- 
action.) However, if only even terms in (a x ) appear in the Hamiltonian, the models have the 
plus-minus symmetry required by Theorem 2.5 of [2] and all of the zeros will lie exactly on the 
unit circle. This shows that the conclusions of the Lee- Yang theorem hold well beyond the set of 
models to which the classic proof applies. 

Finally, in order to demonstrate the non-trivial topology of the set of zeros, let us turn our 
attention to the Blume-Capel model. In this case there are three possible stable phases, each 
corresponding to a particular spin value. In terms of the complex parameter z = e h , the corre- 
sponding metastable free energies are computed from the formulas 

C+(z) = z e A exp {z- 1 e" MJ - A + dz~ 2 ' e -(^2)j-2X + 0(e -4dJ } j ) 

Q_(z) = z-Vexp {ze" MJ - A + dz 2 e-^ d -^ J ' 2X + 0( e - 4dJ )} , (2.25) 

Uz) = exp {(Z + Z-V 2 ^ + d(z 2 + z -2) e -(4d-2)J+2A + ^-^j 

Here it is essential that the energy of the plus-minus neighboring pair exceeds that of zero-plus 
(or zero-minus) by a factor of four. 

A calculation [1] shows that the zeros lie on two curves which are symmetrical with respect to 
circle inversion and which may coincide along an arc of the unit circle, depending on the value 
of A; see Fig. 2. As A increases, the shared portion of these curves grows and, for positive A 
exceeding a constant of order e~ 2dJ , all zeros will lie on the unit circle. Note that by the methods 
of [13], the last result can be established [12] for all temperatures provided A is sufficiently large, 
while our results give the correct critical A but only hold for low temperatures. 

3. Contour models and cluster expansion 

Let us turn to the proofs. We begin by establishing the necessary tools for applying Pirogov-Sinai 
theory. Specifically, we will define contours and show that spin configurations and collections of 
matching contours are in one-to-one correspondence. This will induce a corresponding relation 
between the contour and spin partition functions. We will also summarize the facts we will need 
from the theory of cluster expansions. 

3.1 Contours. 

The goal of this section is to represent spin configurations in terms of contours. Based on the 
fact — following from Assumption C — that the constant configurations are the only possible min- 
ima of (the real part of) the energy, we will define contours as the regions where the spin config- 
uration is not constant. 
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FIGURE 2. A picture demonstrating the location of partition function zeros of the Blume-Capel model. 
Here the zeros concentrate on two curves, related by the circle inversion, which may or may not coincide 
along an arc of the unit circle. There are two critical values of A, denoted by A*, both of order e~ 2dJ , such 
that for A < A^ < 0, the two curves do not intersect; see (a). Once A increases through A<T, a common 
piece starts to develop which grows as A increases through the interval [X~ , A+], see (b) and (c). Finally, 
both curves collapse on the unit circle at A = A+ > and stay there for all A > A+ . With the exception 
of the "bifurcation" points, the zeros lie exactly on the unit circle along the shared arc. The non-uniform 
spacing of the zeros in (b) comes from the influence of the "unstable" phase near the multiple points. 



Recalling our assumption L > 2R + 1, let cr be a spin configuration on IT/, and let Br(<j) 
be the i?-boundary of a. We equip Br(g) with a graph structure by placing an edge between 
any two distinct sites x,y G Br{o~) whenever x and y are contained in a cubic box A C Tl of 
diameter 2R + 1 where a is not constant. We will denote the resulting graph by Gr((t). Some of 
our definitions will involve the connectivity induced by the graph Gr(<j) but we will also use the 
usual concept of connectivity on Tl (or We say that a set of sites A C is connected if 
every two sites from A can be connected by a nearest-neighbor path on A. Note that the connected 
components of Br(<j) and the (vertex sets corresponding to the) components of the graph Gr[ct) 
are often very different sets. 

Now we are ready to define contours. We start with contours on 7L d , and then define contours 
on the torus in such a way that they can be easily embedded into I d . 

Definition 1 A contour on I d is a pair Y = (supply cry) where supp Y is & finite connected 
subset of I d and where cry is a spin configuration on 7h d such that the graph Gr((ty ) is connected 
and Br{o-y) = supp Y. 

A contour on Tl is a pair Y = (supp Y, cry) where supp Y is a non-empty, connected subset 
of Tl with diameter strictly less than L/2 and where cry is a spin configuration on such that 
the graph Gr(cty) is connected and S^(cry) = supp Y. 

A contour network on is a pair N = (supp N, cr^), where N is a (possibly empty or non- 
connected) subset of Tl and where cr^ is a spin configuration on such that Bn{o~ys ) = supp N 
and such that the diameter of the vertex set of each component of Gr(o-jj) is at least L/2. 
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Note that each contour on has an embedding into 7L d which is unique up to translation 
by multiples of L. (Informally, we just need to unwrap the torus without cutting through the 
contour.) As long as we restrict attention only to finite contours, the concept of a contour network 
has no counterpart on Z d , so there we will always assume that N = 0. 

Having defined contours and contour networks onT^ abstractly, our next task is to identify 
the contours Yi, . . . , Y n and the contour network N from a general spin configuration on T^. 
Obviously, the supports of Y\ , . . . , Y n will be defined as the vertex sets of the components of the 
graph Gr{g) with diameter less than L/2, while supp N will be the remaining vertices in Br{g). 
To define the corresponding spin configurations we need to demonstrate that the restriction of a 
to supp Y (resp., suppN) can be extended to spin configurations ay % (resp., er^) on such that 
BR(cTYi) = supp Yi (resp., Br((tj{) = suppN). It will turn out to be sufficient to show that a is 
constant on the boundary of each connected component of Tl \ Br(<j). 

Given a set A C (or A C Z d ), let <9A denote the external boundary of A, i.e., <9A = 
{x G Tl- dist(x, A) = 1}. For the purposes of this section, we also need to define the set A° 
which is just A reduced by the boundary of its complement, A = A \ 8(Tl \ A). An immediate 
consequence of Definition 1 (and the restriction to 2R + 1 > 3) is the following fact: 

Lemma 3.1 Let (A, a) be either a contour or a contour network on T l, and let Cbea connected 
component ofTr. \ A . Then a is constant on C. If (A, a) is a contour on Z d , then a is constant 
on each connected component CofL d \ A , with A now defined as A = A \ d{Z d \ A). 

Proof. Assume that a is not constant on C. Then there must exist a pair of nearest-neighbor 
sites x,y G C such that a x / a y . But then x and all of its nearest neighbors lie in A = Br(g). 
Since C n A = and x G C, we are forced to conclude that x € A \ A°. But that contradicts the 
fact that all of the neighbors of x also lie in A. The same proof applies to contours on Z d . □ 

Definition 2 Let ( A, a) be either a contour or a contour network on Tl and let C be a connected 
component of Tl \ A. The common value of the spin on this component in configuration a 
will be called the label of C. The same definition applies to contours on Z d , and to connected 
components C of Z, d \ A. 

Let A C Tl be a connected set with diameter less than L/2. Since the diameter was defined by 
enclosure into a "cubic" box (see Sect. 1.2), it follows that each such A has a well defined exterior 
and interior. Indeed, any box of side less than L/2 enclosing A contains less than (L/2) d < L d /2 
sites, so we can define the exterior of A, denoted by Ext A, to be the unique component of \ A 
that contains more than L d /2 sites. The interior Int A is defined simply by putting IntA = 
Tl \ (A U Ext A). On the other hand, if A is the union of disjoint connected sets each with 
diameter at least L/2 we define Ext A = and Int A = \ A. These definitions for connected 
sets imply the following definitions for contours on T^: 

Definition 3 Let Y be a contour or a contour network on T^. We then define the exterior 
of Y, denoted by Ext Y, as the set Ext supp Y, and the interior of Y, denoted by Int Y, as the 
set Int supp Y. For each m G S, we let Int m Y be the union of all components of Int Y with 
label m. If Y is a contour on Tl, we say that Y is a m-contour if the label of Ext Y is m. 
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Analogous definitions apply to contours on Z d , except that the exterior of a contour Y is now 
defined as the infinite component of Z d \ supp Y, while the interior is defined as the union of all 
finite components of Z d \ supp Y. 

While most of the following statements can be easily modified to hold for Z d as well as for the 
torus Tl, for the sake of brevity, we henceforth restrict ourselves to the torus. 

Lemma 3.2 Let R > 1 and fix L > 2R + 1. Let a be a spin configuration onTi and let A 
be either the vertex set of a component of the graph Gr{o~) with diameter less than L/2 or the 
union of the vertex sets of all components with diameter at least L/2. Let A' be of the same form 
with A' / A. Then exactly one of the following is true: 

(1) A U Int A C Int A' and A' U Ext A' C Ext A or 

(2) A' U Int A' C Int A and A U Ext A C Ext A', or 

(3) A U Int A C Ext A' and A' U Int A' C Ext A. 

Proof. It is clearly enough to prove the first half of each of the statements (1-3), since the 
second half follow from the first by taking complements (for example in (3), we just use that 
A U Int A C Ext A' implies T L \ (A U Int A) dTl \ Ext A', which is nothing but the statement 
that A' U Int A' C Ext A by our definition of interiors and exteriors). 

In order to prove the first halves of the statements (1-3), we first assume that both A and A' 
are vertex sets of components of the graph Gr{u) with diameter less than L/2. Clearly, since A 
and A' correspond to different components of Gr(<t), we have A n A' = 0. Moreover, A and A' 
are both connected (as subsets of T l) so we have either A C Int A' or A C Ext A' and vice versa. 
Hence, exactly one of the following four statements is true: 

(a) A C Int A' and A' C Int A, or 

(b) A C Int A' and A' C Ext A, or 

(c) A C Ext A' and A' C Int A, or 

(d) A C Ext A' and A' C Ext A. 

We claim that the case (a) cannot happen. Indeed, suppose that A C Int A' and observe that 
if B is a box of size less than L d /2 such that A C B, then Ext A D T L \ B. Hence Int A C B. 
But then B also encloses A and thus Ext AnExt A' D T L \B / 0. Now A'uExt A' is a connected 
set intersecting Ext A but not intersecting A (because we assumed that A C Int A'). It follows 
that A' U Ext A' C Ext A, and hence Int A'dAU Int A. But then we cannot have A' C Int A as 
well. This excludes the case (a) above, and also shows that (b) actually gives A U Int A C Int A', 
which is the first part of the claim (1), while (c) gives A' U Int A' C Int A, which is the first part 
of the claim (2). 

Turning to the remaining case (d), let us observe that A' C Ext A implies Int An A' C 
Int A n Ext A = 0. Since A n A' = as well, this implies (A U Int A) n A' = 0. But 
A U Int A is a connected subset of T L , so either A U Int A C Int A' or A U Int A C Ext A'. 
Since A C Ext A' excludes the first possibility, we have shown that in case (d), we necessarily 
have A U Int A C Ext A', which is the first part of statement (3). This concludes the proof of the 
lemma for the case when both A and A' are vertex sets of components of the graph Gn(a) with 
diameter less than L/2. 
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Since it is not possible that both A and A' are the union of the vertex sets of all components 
of diameter at least L/2, it remains to show the statement of the lemma for the case when A is 
the vertex set of a component of the graph Gr{o~) with diameter less than L/2, while A' is the 
union of the vertex sets of all components of diameter at least L/2. By definition we now have 
Ext A' = 0, so we will have to prove that A U Int A C Int A', or equivalently, A' C Ext A. 
To this end, let us first observe that A n A' = 0, since A has diameter less than L/2 while all 
components of A' have diameter at least L/2. Consider the set Int A. Since A has diameter less 
than L/2, we can find a box B of side length smaller than L/2 that contains A, and hence also 
Int A. But this implies that none of the components of A' can lie in Int A (their diameter is too 
large). Since all these components are connected subsets of Int A U Ext A, we conclude that they 
must be part of Ext A. This gives the desired conclusion A' C Ext A. □ 

The previous lemma allows us to organize the components of Gr(ct) into a tree-like structure 
by regarding A' to be the "ancestor" of A (or, equivalently, A to be a "descendant" of A') if the 
first option in Lemma 3.2 occurs. Explicitly, let Wr(g) be the collection of all sets A C Tl 
that are either the vertex set of a connected component of Gr{o~) with diameter less than L/2 or 
the union of the vertex sets of all connected components of diameter at least L/2. We use Ao to 
denote the latter. If there is no component of diameter L/2 or larger, we define Ao = and set 
IntA = T L . 

We now define a partial order on Wr(g) by setting A -< A' whenever A U Int A C Int A'. 
If A -< A', but there is no A" G W R (a) such that A -< A" -< A', we say that A is a child 
of A' and A' is a parent of A. Using Lemma 3.2, one easily shows that no child has more than 
one parent, implying that the parent child relationship leads to a tree structure on Wr{o~), with 
root Ao- This opens the possibility for inductive arguments from the innermost contours (the 
leaves in the above tree) to the outermost contours (the children of the root). Our first use of such 
an argument will be to prove that unique labels can be assigned to the connected components of 
the complement of Br{o~) . 

Lemma 3.3 Let a be a spin configuration onT^ and let A be either the vertex set of a component 
of the graph Gr(ct) with diameter less than L/2 or the set of sites in Br(g) that are not contained 
in any such component. If C is a connected component ofTr, \ A , then a is constant on C PI A. 

The proof is based on the following fact which is presumably well known: 

Lemma 3.4 Let A C Z rf be a finite connected set with a connected complement. Then dA c is 
^-connected in the sense that any two sites x, y £ dA c are connected by a path on dA c whose 
individual steps connect only pairs of sites of7L d with Euclidean distance not exceeding \[2. 

Proof. The proof will proceed in three steps. In the first step, we will prove that the edge boundary 
of A, henceforth denoted by 5 A, is a minimal cutset. (Here we recall that a set of edges E' in a 
graph G = (V, E) is called a cutset if the graph G' = (V,E\ E') has at least two components, 
and a cutset E' is called minimal if any proper subset of E' is not a cutset.) In the second step, 
we will prove that the dual of the edge boundary 5 A is a connected set of facets, and in the third 
step we will use this fact to prove that dA c is *-connected. 

Consider thus a set A which is connected and whose complement is connected. Let 5 A be the 
edge boundary of A and let Ed be the set of nearest-neighbor edges in Z d . The set 5 A is clearly a 
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cutset since any nearest-neighbor path joining A to A c must pass through one of the edges in 5 A. 
To show that 5 A is also minimal, let E' be a proper subset of 5 A, and let e G 5 A \ E' . Since 
both A and A c are connected, an arbitrary pair of sites x, y G Z, d can be joined by a path that uses 
only edges in {e} U (Ed \ 5 A) C \ E'. Hence such E' is not a cutset which implies that 8 A 
is minimal as claimed. 

To continue with the second part of the proof, we need to introduce some notation. As usual, 
we use the symbol Z* d to denote the set of all points in M. d with half-integer coordinates. We 
say that a set c C Z* d is a k-cell if the vertices in c are the "corners" of a /c-dimensional unit 
cube in R d . A d-cell c C Z* d and a vertex x G Z d are called dual to each other if x is the 
center of c (considered as a subset of R d ). Similarly, a facet / (i.e., a (d — l)-cell in Z* d ) and a 
nearest-neighbor edge e C Z rf are called dual to each other if the midpoint of e (considered as a 
line segment in M. d ) is the center of /. The boundary dC of a set C of d-cells in Z* d is defined 
as the set of facets that are contained in an odd number of cells in C, and the boundary dF of a 
set F of facets in Z* d is defined as the set of (d — 2)-cells that are contained in an odd number of 
facets in F. Finally, a set of facets F is called connected if any two facets /, /' G -F can be joined 
by a path of facets /i = /,...,/„ = /' in F such that for all z = 1, . . . , n — 1, the facets fi 
and f i+ i share a (d - 2)-cell in Z* d . 

Note that an arbitrary finite set of facets F has empty boundary if and only if there exists a 
finite set of cubes C such that F = dC, which follows immediately from the fact R d has trivial 
homology. Using this fact, we now prove that the set F of facets dual to 5 A is connected. Let 
W be the set of d-cells dual to A, and let F = dW be the boundary of W . We will now prove 
that F is a connected set of facets. Indeed, since F = dW, we have that F has empty boundary, 
dF = 0. Assume that F has more than one component, and let F C F be one of them. Then F 
and F \ F are not connected to each other, and hence share no (d — 2)-cells. But this implies that 
the boundary of F must be empty itself, so that F is the boundary of some set W. This in turn 
implies that the dual of F is a cutset, contradicting the fact that 5 A is a minimal cutset. 

Consider now two points x, y G dA° C A. Then there are points x, y G A° such that {x, x} 
and {y, y} are edges in 5 A. Taking into account the connectedness of the dual of 5 A, we can find 
a sequence of edges e\ = {x, x}, e n = {y, y} in 5 A such that for all k = 1, . . . , n — 1, the 
facets dual to and e^+i share a (d— 2) cell in Z* d . As a consequence, the edges and eu+i are 
either parallel, and the four vertices in these two edges form an elementary plaquette of the form 
{x, x+ni, x+ri2, x+ni+112} where ni and ni are unit vectors in two different lattice directions, 
or efc and e^+i are orthogonal and share exactly one endpoint. Since both and e£. +1 are edges 
in 5 A, each of them must contains a point in 8A°, and by the above case analysis, the two points 
are at most y/2 apart. The sequence ei, . . . , e n thus gives rise to a sequence of (not necessarily 
distinct) points xi, . . . , x n G dA c such that x = x\, y = x n and dist(xfc, Xk+i) < for all 
k = 1, . . . ,n — 1. This proves that dA c is * -connected. □ 

Proof of Lemma 3.3. Relying on Lemma 3.2, we will prove the statement by induction from 
innermost to outermost components of diameter less than L/2. Let A be the vertex set of a 
component of the graph Gr(<t) with diameter less than L/2 and suppose Bji(a) n Int A = 0. (In 
other words, A is an innermost component of Bn(a).) Then the same argument that was used in 
the proof of Lemma 3.1 shows that all connected components of Int A clearly have the desired 
property, so we only need to focus on Ext A. 
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Let us pick two sites x,y G d Ext A = An d Ext A and let A' = A U Int A. Then A' is 
connected with a connected complement and since A has a diameter less than L/2, we may as 
well think of A' as a subset of Z rf . Now Lemma 3.4 guarantees that <9(A') C = d Ext A is *- 
connected and hence x and y are connected by a *-connected path entirely contained in d Ext A. 
But the spin configuration must be constant on any box (z + [-R, R] d ) fl Z d with z G d Ext A 
and thus the spin is constant along the path. It follows that a x = a y . 

The outcome of the previous argument is that now we can "rewrite" the configuration on A' 
without changing the rest of Br(<t). The resulting configuration will have fewer connected com- 
ponents of diameter less than L/2 and, proceeding by induction, the proof is reduced to the cases 
when there are no such components at all. But then we are down to the case when A simply 
equals Bji(a). Using again the argument in the proof of Lemma 3.1, the spin must be constant 
on each connected component C of \ Br(o-)°. □ 

The previous lemma shows that each component of the graph Gji(a) induces a unique label 
on every connected component C of its complement. Consequently, if two contours share such a 
component — which includes the case when their supports are adjacent to each other — they must 
induce the same label on it. A precise statement of this "matching" condition is as follows. (Note, 
however, that not all collections of contours will have this matching property.) 

Definition 4 We say that the pair (Y, N) — where Y is a set of contours and N is a contour 
network on — is a collection of matching contours if the following is true: 

(1) suppY n suppY' = for any two distinct Y, Y' G Y and suppY n suppN" = for 
any Y G Y. 

(2) If C is a connected component of \ [(suppN) U Uy e Y( su PP^)°]' tnen tne restrictions 
of the spin configurations ay (and <r^) to C are the same for all contours Y G Y (and contour 
network K) with supp Y n C / (supp N fl C ^ 0). In other words, the contours/contour 
network intersecting C induce the same label on C. 

Here we use the convention that there are altogether |<S| distinct pairs (Y, K) with both Y = 
and N = 0, each of which corresponds to one m G S. 

Definition 4 has an obvious analogue for sets Y of contours on 7j d , where we require that 
(1) suppY n suppY' = for any two distinct Y,Y' G Y and (2) all contours intersecting a 
connected component C of 7L d \ [Uy£y( su PP^)°] induce the same label on C. 

It remains to check the intuitively obvious fact that spin configurations and collections of 
matching contours are in one-to-one correspondence: 

Lemma 3.5 For each spin configuration a G S Th , there exists a unique collection (Y, N) of 
matching contours onT^ and for any collection (Y, N) of matching contours on T^, there exists 
a unique spin configuration a G <S Tl such that the following is true: 

(1) The supports of the contours in Y (of the contour network NJ are the vertex sets (the union of 
the vertex sets) of the connected components of the graph Gn(a) with diameter strictly less 
than (at least) L/2. 

(2) The spin configuration corresponding to a collection (Y, K) of matching contours arise by 
restricting oy for each Y G Y as well as to the support of the corresponding contour 
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(contour network) and then extending the resulting configuration by the common label of the 
adjacent connected components. 

Proof. Let a be a spin configuration and let A be a component of the graph Gr((t) with diameter 
less than L/2. Then Lemma 3.3 ensures that a is constant on the boundary dC of each component 
C of A c . Restricting a to A and extending the resulting configuration in such a way that the new 
configuration, a, restricted to a component component C of A c , is equal to the old configuration 
on dC, the pair (A, a) thus defines a contour. Similarly, if A is the union of all components of 
the graph Gr{o~) with diameter at least L/2 and C is a connected component of \ A°, then a 
is, after removal of all contours, constant on C. The contours/contour network (Y, K) then arise 
from a in the way described. The supports of these objects are all disjoint, so the last property to 
check is that the labels induced on the adjacent connected components indeed match. But this is 
a direct consequence of the construction. 

To prove the converse, let (Y, N) denote a set of matching contours and let a be defined by 
the corresponding contour configuration on the support of the contours (or contour network) and 
by the common value of the spin in contour configurations for contours adjacent to a connected 
component of \ [(suppN) U |Jy eY (supp Y)°]. (If at least one of Y, IN" is nonempty, then 
this value is uniquely specified because of the matching condition; otherwise, it follows by our 
convention that empty (Y, N) carries an extra label.) 

It remains to show that Y are the contours and N is the contour network of a. Let A be 
a component of the graph Gr(ct). We have to show that it coincides with suppY for some 
Y G Y or with a component of suppDsf (viewed as a graph). We start with the observation 
that A C supp IN" U (UreY supp Y). Next we note that for each Y G Y, the graph Gu(ay) is 
connected. Since the restriction of ay to supp Y is equal to the corresponding restriction of a, 
we conclude that supp Y n A / implies supp Y C A, and similarly for the components of 
supp INT. To complete the proof, we therefore only have to exclude that supp Y C A for more 
than one contour Y G Y, or that A C A for more than one component A of supp N, and similarly 
for the combination of contours in Y and components of supp N. 

Let us assume that supp Y C A for more than one contour Y G Y. Since A is a connected 
component of the graph Gr{g), this implies that there exists a box B z = (z + [-R, R] d ) H T, d 
and two contours Y±, Y2 € Y such that a is not constant on B z , supp Yi U supp Y2 C A and B z 
is intersecting both suppYi and suppY2. But this is in contradiction with the fact that Y is 
a collection of matching contours (and a configuration on any such box not contained in the 
support of one of the contours in Y or in a component of supp IN" must be constant). In the same 
way one excludes the case combining supp Y with a component of supp IN or combining two 
components of supp IN. Having excluded everything else, we thus have shown that A is either 
the support of one of the contours in Y, or one of the components of supp IN. □ 

3.2 Partition functions and Peierls' condition. 

A crucial part of our forthcoming derivations concerns various contour partition functions, so our 
next task will be to define these quantities. We need some notation: Let (Y, N) be a collection 
of matching contours on Tl. A contour Y G Y is called an external contour in Y if supp Y C 
Ext Y 1 for all Y' G Y different from Y, and we will call two contours Y,Y' G Y mutually 
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external if supp Y C Ext Y' and supp Y' C Ext Y. Completely analogous definitions apply 
to a set of matching contours Y on 7L d (recall that on Z d , we always set N = 0). Note that, 
by Lemma 3.2, two contours of a configuration a on T^ are either mutually external or one is 
contained in the interior of the other. Inspecting the proof of this Lemma 3.2, the reader may 
easily verify that this remains true for configurations on Z d , provided the set Br(g) is finite. 

Given a contour Y = (supp Y, ay) or a contour network N = (supp N, cr^) let E(Y,z) 
and E(N, z) denote the corresponding excitation energies E(ay, z) and E(aj^, z) from (1.7). We 
then introduce exponential weights p z (Y) and p z (N), which are related to the quantities E(Y, z) 
and E(N, z) according to 

p z (Y) = e~ E ^ and p z (Ji) = e~ E ^' z \ (3.1) 

The next lemma states that the exponential weights d m (z), p z {Y) an d Pz(N) are analytic func- 
tions of z. 

Lemma 3.6 Suppose that Assumption CO holds, let q G S, let Y be a q-contour and let N be a 
contour network. Then 6 q (z), p z (Y) and p z (N) are analytic functions of z in G. 

Proof. By assumption CO, the functions z ^ ip\(a,z) = exp{— $a(ct, z)} are holomorphic 
in &. To prove the lemma, we will show that 9 q (z), p z (Y) and p z (N) can be written as products 
over the exponential potentials ip\{o-, z), with a = a q , a = ay and a = a^, respectively. 

Let us start with 9 q (z). Showing that 9 q is the product of exponential potentials (fA(a q , z) is 
clearly equivalent to showing that e q can be rewritten in the form 

e q = $ aK), (3-2) 

Aev e 

where V e is a collection of subsets A C Tl- But this is obvious from the definition (1.6) of e q : 
just choose V e in such a way that it contains exactly one representative from each equivalence 
class under translations. 

Consider now a contour Y = (supp Y, cry) and the corresponding excitation energy E(Y, z). 
We will want to show that E(Y, z) can be written in the form 

E(Y,z) = $aM, (3.3) 
AeVy 

where Vy is again a collection of subsets A C Tl. Let A q = Ext Y U lnt q Y, and A m = Int m Y 
for m q. Consider a point x € A m . Since x ^ suppy = Br{o~y), the configuration ay must 
be constant on any subset A C Tl that has diameter 2R + 1 or less and contains the point x, 
implying that 

E m^)= E m**{° m ) = em (3-4) 

A: xGA 11 A: xGA 1 1 
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whenever x G A m . Using these facts, we now rewrite E(Y, z) as 

e(y, z) = 0H L (a Y ) EE wM a ^ 
= E - E i A ™i e ™ 

AcT L meS (3-5) 

= £ 

AcsuppY meS I v AcT L 

AnA m ^0 

To complete the proof, we note that the sum over all A with A n A m / contains at least |A m | 
translates of each A C T/, contributing to the right hand side of (3.2). As a consequence, the 
difference on the right hand side of (3.5) can be written in the form (3.3), proving that E(Y, z) is 
of the form (3.3). The proof that p z (J{) is an analytic function of z is virtually identical. □ 

Next we define partition functions in finite subsets of Z d . Fix an index q G S. Let A C T, d 
be a finite set and let M(A, q) be the set of all collections Y of matching contours in 7L d with the 
following properties: 

(1) For each Y G Y, we have supp Y U Int Y C A. 

(2) The external contours in Y are g-contours. 

Note that supp Y U Int Y C A is implied by the simpler condition that supp Y C A if Z d \ A is 
connected, while in the case where Z, d \ A is not connected, the condition supp Y U Int Y C A 
is stronger, since it implies that none of the contours Y G Y contain any hole of A in its interior. 
(Here a hole is defined as a finite component of Z d \ A.) In the sequel, we will say that Y is a 
contour in A whenever Y obeys the condition supp Y U Int Y C A. 

The contour partition function in A with boundary condition q is then defined by 

z q (\, Z )= [ n M^) |Am(¥)i ] n ^ 

YeM(A,q) meS YGY 

where A m (Y) denotes the union of all components of A \ Uygy su PP^ w ^ tn m, and 
|A m (Y)| stands for the cardinality of A m (Y). 

If we add the condition that the contour network N is empty, the definitions of the set Ai (A, q) 
and the partition function Z q (A, z) clearly extends to any subset A C Tl, because on every 
contour has a well defined exterior and interior. However, our goal is to have a contour repre- 
sentation for the full torus partition function. Let M l denote the set of all collections (Y, Ji) of 
matching contours in Tl which, according to our convention, include an extra label meS when 
both Y and N are empty. If (Y, N) G Ml is such a collection, let A m (Y, N) denote the union of 
the components of Tl \ (supp N U Uygy su PP Y) with label m. Then we have: 

Proposition 3.7 (Contour representation) The partition function on the torus T l is given by 

z t^)= e [n^) |AmOT ]fe(N)nfe(^)' (3-7) 
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In particular, we have 

ZT v {z)= Y, P'W II Z m (A m ((H,J{),z). (3.8) 
(8,N)eM L mes 

Proof. By Lemma 3.5, the spin configurations a are in one-to-one correspondence with the pairs 
(Y, N) G Ml- Let (Y, N) be the pair corresponding to a. Rewriting (1.8) as 

f3H L {a) = Y E m $A ^)' (3 " 9) 

seT L A: ACT L 1 1 
A3x 

we can now split the first sum into several parts: one for each m G S corresponding to x G 
A m (Y,N), one for each Y G Y corresponding to a; G suppF, and finally, one for the part of 
the sum corresponding to x G suppN. Invoking the definitions of the energies e m (z), E(Y, z) 
and E(N, z), this gives 

PH L (a)= Y e m(z)\A m (Y,N)\ + Y E ( Y ' z )+ E (K,z). (3.10) 

meS YeY 

Strictly speaking, the fact that the excitation energy factors (technically, sums) over contours and 
contour networks requires a proof. Since this is straightforward using induction as in the proof of 
Lemma 3.3, starting again with the innermost contours, we leave the formal proof to the reader. 
Using the definitions of 6 m (z), p z (Y) and p z (N) and noting that, by Lemma 3.5, the sum over a 
can be rewritten as the sum over (Y, N) G M-l, formula (3.7) directly follows. 

The second formula, (3.8), formally arises by a resummation of all contours that can contribute 
together with a given contour network N. It only remains to check that if Y m C Y is the set 
of Y G Y with supp Y C A m = A m (0, X), then Y m can take any value in M (A m , m). But this 
follows directly from Definition 4 and the definition of M (A m , m). □ 

In order to be useful, the representations (3.7) and (3.8) require that contours and contour 
networks are sufficiently suppressed with respect to the maximal ground state weight 6. This 
is ensured by Assumption C2, which guarantees that |p z (y)| < 6>(z)! y 'e~ r l y l and < 
9(z)^e~ T ^\, where we used the symbols \Y\ and |N| to denote the cardinality of suppF and 
supp N, respectively. 

3.3 Cluster expansion. 

The last ingredient that we will need is the cluster expansion, which will serve as our principal 
tool for evaluating and estimating logarithms of various partition functions. The cluster expansion 
is conveniently formulated in the context of so-called abstract polymer models [19, 10, 7, 14]. 
Let K be a countable set — the set of all polymers — and let ^ be the relation of incompatibility 
which is a reflexive and symmetric binary relation on K. For each A C K, let M(A) be the set 
of multi-indices X: K — > {0} U N that are finite, ^ 76K X(7) < oo, and that satisfy X(7) = 
whenever 7 A. Further, let C(A) be the set of all multi-indices X G M. (A) with values in {0, 1} 
that satisfy X(7)X(7') = whenever 7 9^ 7' and 7 / 7'. 
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Let 3 : K — > C be a polymer functional. For each finite subset A C K, let us define the polymer 
partition function 2 (A) by the formula 

Z (A)= E I[3(7) X(7) - (3-11) 

XGC(A) 7GK 

In the most recent formulation [7, 14], the cluster expansion corresponds to a multidimensional 
Taylor series for the quantity log -2(A), where the complex variables are the 3(7). Here clusters 
are simply multi-indices X € M(K) for which any nontrivial decomposition of X leads to incom- 
patible multi-indices. Explicitly, if X can be written as Xi + X2 with Xi, X2 ^ 0, then there exist 
two (not necessary distinct) polymers 71, 72 £ K, 71 9^ 72, such that Xi (7^X2(72) 7^ 0. 

Given a finite sequence r = (71, . . . , j n ) of polymers in K, let n(T) = n be the length of the 
sequence Y, let Q{T) be the set of all connected graphs on {1, . . . , n} that have no edge between 
the vertices i and j if 7^ ~ jj, and let Xr be the multi-index for which Xr(7) is equal to the 
number of times that 7 appears in V. For a finite multi-index X, we then define 

« T w= E ^ E ( 3 - 12 ) 

r:x r =x v ; seg(r) 
with \g\ denoting the number of edges in Q, and 

3 T (X) = « T (X)[l3(7) X(7) . (3-13) 

7GK 

Note that Q(T) = if Xr is not a cluster, implying, in particular, that 3 T (X) = whenever X is 
not a cluster. We also use the notation X 9^ 7 whenever X is a cluster such that X(7 ; ) > for at 
least one 7' 9^ 7. 

The main result of [14] (building upon [7]) is then as follows: 

Theorem 3.8 (Cluster expansion) Let a : K — > [0, 00) be a function and let 30 : K — > [0, 00) be 
polymer weights satisfying the bound 

^3o(7V (y) <a( 7 ), 7GK. (3.14) 

7'€K 

r/ie« 2(A) 7^ 0/or any finite set A C K a«J a«j collection of polymer weights 3 : K — > C in the 
multidisc D A = {(3(7)) : |3( 7 )| < 3o(7)i 7 £ A}. Moreover, if we define log 2(A) as unique 
continuous branch of the complex logarithm of 2(A) o« D A normalized so that log 2(A) = 
w/ie« 3(7) = Ofor all 7 G A, ?/ie« 

log 2(A) = J] 3 T (X) (3.15) 
xeX(A) 

holds for each finite set A C K. //ere ?/ie power series on the right hand side converges absolutely 
on the multidisc D A . Furthermore, the bounds 

E |3 T (X)| < E X( 7 )|3 T (X)| < |3(7)|e a(7) (3.16) 
xe.M(K) xeM(K) 

X(7)>1 
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and 

|3 T (X)|<a( 7 ) (3.17) 

xeM(K) 

X^ 7 

hold for each 7 G K. 

Proof. This is essentially the main result of [14] stated under the (stronger) condition (3.14), 
which is originally due to [15, 10]. To make the correspondence with [14] explicit, let 

A*( 7 ) = log(l + |3( 7 )|e a W) (3.18) 

and note that ^(7) < |3(7)|e a ( 7 ) < $o(j)e a (~<\ The condition (3.14) then guarantees that we 
have Sy^MY) < 0(7) and hence 

13(7)1 = (e M(7) - l)e" a(7) < (e"W - l)exp{- £ M ( 7 ')}. (3.19) 

7' 7^7 

This implies that any collection of weights 3 : K — ► C such that (3(7)1 < 3o( 7 ) f° r all 7 G K will 
fulfill the principal condition of the main theorem of [14]. Hence, we can conclude that Z(A) 7^ 
in D A and that (3.15) holds. Moreover, as shown in [14], both quantities on the left-hand side of 
(3.16) are bounded by — 1 which simply equals |3 (t) I e -^^ . The bound (3.16) together with 
the condition (3.14) immediately give (3.17). □ 

To facilitate the future use of this result, we will extract the relevant conclusions into two 
lemmas. Given a spin state q G S, let K q denote the set of all g-contours in Z d . If Y, Y' G K q , let 
us call Y and Y 1 incompatible if supp Y n supp Y' / 0. If A is a finite set of (/-contours, we will 
let Z(A) be the polymer sum (3.11) defined using this incompatibility relation. Then we have: 

Lemma 3.9 There exists a constant cq = co(d, |<S|) G (0, 00) such that, for all q G S and all 
contour functionals 3 : K q — ► C satisfying the condition 

\i{Y)\ < 3o0n = e-( C0 +")l y l for all Y G K q , (3.20) 
for some rj > 0, the following holds for all k > 1: 

(1) Z(A) / for all finite Ac/C, with log Z(A) given by (3.15), and 

Y, |3 T (X)|<e^ fc . (3.21) 

X£M(Kq) 

V(X)30, ||X||>fc 

HereV{X) = \J Y . x(y)>0 V(Y) withV{Y) = suppF U IntY and ||X|| = Ey eK9 X ( Y )\ Y \- 

(2) Furthermore, if the activities %(Y) are twice continuously differ entiable (but not necessarily 
analytic) functions of a complex parameter z such that the bounds 

\d w l{Y)\ <3o(y) and \d w d w ,}(Y)\ < io (Y) (3.22) 

hold for any w, w' G {z, z} and any Y G K q , then 

Y l 5 -3 T (X)| <e- r ' k and |^^3 T (X) | < e"^. (3.23) 

X£M(K q ) X£M(K q ) 
V(X)30,||X||>fc V{X)B0, \\X\\>k 
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for any w, w' G {z, z}. 

Using, for any finite A C Z d , the notation K 9; a = {Y G K q : supp Y U Int Y C A} and <9A 
for the set of sites in Z d \ A that have a nearest neighbor in A, we get the following lemma as an 
easy corollary: 

Lemma 3.10 Suppose that the weights 3 satisfy the bound (3.20) and are invariant under the 
translations ofZ d . Then the polymer pressure s q = lim A | Z d | A| _1 log Z(K qt A) exists and is 
given by 

s q = E n7^ T (X). (3.24) 



xe.vi(K„): \ ixeo ^ X )' 



Moreover, the bounds 
and 



s q \ < e~ n (3.25) 



|logZ(K,, A ) - »,|A|| < e-"|8A| (3.26) 
hold. Finally, if the conditions (3.22) on derivatives of the weights ^(Y) are also met, the polymer 
pressure s q is twice continuously differentiable in z with the bounds 

\d w s q \ < e~' ] and \d w d w is q \ < e~ v , (3.27) 

valid for any w, w' £ {z, z}. 

Proof of Lemma 3. 9. Let us consider a polymer model where polymers are either a single site 
of Z d or a (/-contour from K^. (The reason for including single sites as polymers will become 
apparent below.) Let the compatibility between contours be defined by disjointness of their sup- 
ports while that between a contour Y and a site x by disjointness of {x} and supp Y U Int Y. If 
we let 3(7) = whenever 7 is just a single site, this polymer model is indistinguishable from the 
one considered in the statement of the lemma. Let us choose cq so that 

e ( 2 - c °)l y l < 1. (3.28) 

YeK q : V(Y)3Q 

To see that this is possible with a constant c depending only on the dimension and the cardinality 
of S, we note that each polymer is a connected subset of Z d . As is well known, the number of 
such sets of size n containing the origin grows only exponentially with n. Since there are only 
finitely many spin states, this shows that it is possible to choose cq as claimed. 

Defining 0(7) = 1 if 7 is a single site and a(Y) = \Y\ if Y is a g-contour in K q , the assump- 
tion (3.14) of Theorem 3.8 is then satisfied. (Note that assumption (3.14) requires slightly less 
than (3.28), namely the analogue of (3.28) with the exponent of (1 — cn)|Y| instead of (2 — cn)|Y|; 
the reason why we chose cq such that (3.28) holds will become clear momentarily.) Theorem 3.8 
guarantees that Z (A) 7^ and (3.15) holds for the corresponding cluster weights 3 T . Actually, 
assumption (3.14) is, for all 77 > 0, also satisfied when j(Y) is replaced by $(Y)e b ( Y ^ with 
b(Y) = rj\Y\, yielding 

E e 6(X) |3 T (X)|<a( 7 ) (3.29) 
xeA4(K) 

X^7 
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with b(X) = rj\\X\\ instead of (3.17). Using (3.29) with 7 chosen to be the polymer represented by 
the site at the origin and observing that the quantity 6(X) exceeds r]k for any cluster contributing 
to the sum in (3.21), we get the bound 

c"* E |3 T (X)|< £ |3 T (X)|e 6 W < 1, (3.30) 
v(x)30, ||x||>fc y(x)90 

i.e., the bound (3.21). 

In order to prove the bounds (3.23), we first notice that, in view of (3.13) and (3.22) we have 

\d w f(X) I < ||X|| | 3 5(X) I < e" x ll | 3 J(X) I (3.31) 

and 

\d w d w ,f(X)\ < ||X|| 2 |3j(X)| < e i|x|l |jJ(X)|. (3.32) 

Using (3.29) with b(Y) = (77 + 1) | V| (which is also possible since we choose Co such that (3.28) 
holds as stated, instead of the weaker condition where (2 — co)|^| is replaced by (1 — co)|^|) we 
get (3.23) in the same way as (3.21). □ 

Proof of Lemma 3.10. The bound (3.21) for k = 1 immediately implies that the sum in (3.24) 
converges with \s q \ < e~ v . Using (3.15) and standard resummation techniques, we rewrite the 
left hand side of (3.26) as 



|logZ ( K g>A )- Sg |A|| = | £ ^^V(X) 

v(x)£A 



(3.33) 



Next we note that for any cluster X G M(K q ), the set V(X) is a connected subset of Z d , which 
follows immediately from the observations that supp Y U Int Y is connected for all contours Y, 
and that incompatibility of two contours Y, Y' implies that supp Y n supp Y' 7^ 0. Since only 
clusters with V(X) n A / and V(X) n A c / contribute to the right hand side of (3.33), we 
conclude that the right hand side of (3.33) can be bounded by a sum over clusters X G K q with 
V(X) n dA / 0. Using this fact and the bound (3.21) with k = 1, (3.26) is proved. 

Similarly, using the bounds (3.23) in combination with explicit expression (3.24) in terms of 
absolutely converging cluster expansions, the claims (3.27) immediately follow. □ 

Remark 5. The proof of Lemma 3.9 holds without changes if we replace the set of all ^-contours 
in Z d by the set of all ^-contours on the torus Tr,. This is not true, however, for the proof of the 
bound (3.26) from Lemma 3.10 since one also has to take into account the difference between 
clusters wrapped around the torus and clusters in Z d . The corresponding modifications will be 
discussed in Section 4.4. 



4. PlROGOV-SlNAI ANALYSIS 



The main goal of this section is to develop the techniques needed to control the torus partition 
function. Along the way we will establish some basic properties of the metastable free energies 
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which will be used to prove the statements concerning the quantities ( m . Most of this section 
concerns the contour model on Z d . We will return to the torus in Sections 4.4 and 5. 

All of the derivations in this section are based on assumptions that are slightly more general 
than Assumption C. Specifically, we only make statements concerning a contour model satisfying 
the following three conditions (which depend on two parameters, r and M): 

(1) The partition functions Z q (A,z) and Z£ cr (z) are expressed in terms of the energy vari- 
ables 9 m (z) and contour weights p z as stated in (3.6) and (3.7), respectively. 

(2) The weights p z of contours and contour networks are translation invariant and are twice 
continuously differentiable functions on G. They obey the bounds 

\dldlp z (Y)\ < (M\Y\) e+i e-^9(z)^, (4.1) 

and 

\didl Pz (Ji)\ < (M\X\y +l e-^e(z)W (4.2) 

as long as £, £ > and £ + £ < 2. 

(3) The energy variables 9 m wee twice continuously differentiable functions on G and obey the 
bounds 

\d e z dl9 m (z)\ < (M) e+i 6(z) (4.3) 

as long as £, I > and £ + £ < 2. We will assume that 9(z) is bounded uniformly from below 
throughout G. However, we allow that some of the 9 m vanish at some z G G. 

In particular, throughout this section we will not require that any of the quantities 9 m , p z (Y) 
or |0 2 (N) is analytic in z. 

4.1 Truncated contour weights. 

The key idea of contour expansions is that, for phases that are thermodynamically stable, contours 
appear as heavily suppressed perturbations of the corresponding ground states. At the points of 
the phase diagram where all ground states lead to stable phases, cluster expansion should then 
allow us to calculate all important physical quantities. However, even in these special circum- 
stances, the direct use of the cluster expansion on (3.6) is impeded by the presence of the energy 
terms 9 m {z)\ K ™^\ and, more seriously, by the requirement that the contour labels match. 

To solve these problems, we will express the partition function in a form which does not involve 
any matching condition. First we will rewrite the sum in (3.6) as a sum over mutually external 
contours Y ext times a sum over collections of contours which are contained in the interior of one 
of the contours in Y ext . For a fixed contour Y € Y ext , the sum over all contours inside Int m Y 
then contributes the factor Z m (lnt m Y, z), while the exterior of the contours in Y ext contributes 
the factor m (z)\ Ext l, where Ext = Ext A (Y ext ) = n YeY =x.(Ext Y n A). As a consequence, we 
can rewrite the partition function (3.6) as 

Z q (A,z) = ^9 q (z)\ Ext \ H [ Pz (Y)l[Z m (Int m Y,z)}, (4.4) 
where the sum goes over all collections of compatible external ^-contours in A. 
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At this point, we use an idea which originally goes back to [9]. Let us multiply each term in 
the above sum by 1 in the form 

Z q (lnt m Y,z) 



n n 



Z a (Int m Y, z) 



(4.5) 



Associating the partition functions in the denominator with the corresponding contour, we get 

Z q {K,z) = Y,O q {z)\^\ [J (0 q (z)\ Y \K q (Y,z)Z q (lntY,z)), (4.6) 

where K q (Y, z) is given by 

K q (Y, z) = Pz (Y) O q (z)-W J] IWW- ^ 
Proceeding by induction, this leads to the representation 

Z q (A,z)=e q (z)\ A \ £ H K <i( Y > z )> ^ 

YeC(A,q) YeY 

where C(A, q) denotes the set of all collections of non-overlapping (/-contours in A. Clearly, the 
sum on the right hand side is exactly of the form needed to apply cluster expansion, provided the 
contour weights satisfy the necessary convergence assumptions. 

Notwithstanding the appeal of the previous construction, a bit of caution is necessary. Indeed, 
in order for the weights K q (Y, z) to be well defined, we are forced to assume — or prove by 
cluster expansion, provided we somehow know that the weights K q have the required decay — that 
Z q (lnt m Y, z) / 0. In the "physical" cases when the contour weights are real and positive (and 
the ground-state energies are real-valued), this condition usually follows automatically. However, 
here we are considering contour models with general complex weights and, in fact, our ultimate 
goal is actually to look at situations where a partition function vanishes. 

Matters get even more complicated whenever there is a ground state which fails to yield a stable 
state (which is what happens at a generic point of the phase diagram). Indeed, for such ground 
states, the occurrence of a large contour provides a mechanism for flipping from an unstable to 
a stable phase — which is the favored situation once the volume gain of free energy exceeds the 
energy penalty at the boundary. Consequently, the relative weights of (large) contours in unstable 
phases are generally large, which precludes the use of the cluster expansion altogether. A classic 
solution to this difficulty is to modify the contour functionals for unstable phases [21, 5, 6]. We 
will follow the strategy of [6], where contour weights are truncated with the aid of a smooth 
mollifier. 

To introduce the truncated contours weights, let us consider a C 2 (R) -function x i-> xi x )> sucn 
that < %(■) < 1, x( x ) = for x < —2 and x{ x ) = 1 for x > —1. Let cq be the constant from 
Lemma 3.9. Using x as a regularized truncation factor, we shall inductively define new contour 
weights K' q (-,z) so that \K' q (Y,z)\ < e"( C0+r / 2 )l y l for all g-contours Y. By Lemma 3.9, the 
associated partition functions Z' q (-,z) defined by 

Z>(A,z)=9 q (z)\ A \ £ l[K> q (Y,z) (4.9) 

YeC(A, 9 ) yeY 



PARTITION FUNCTION ZEROS AT FIRST-ORDER PHASE TRANSITIONS 



29 



can then be controlled by cluster expansion. (Of course, later we will show that K' q (-,z) = 
K q (-, z) and Z' q (A, z) = Z q (A, z) whenever the ground state q gives rise to a stable phase.) 

Let 9 q (z) 7^ 0, let Y be a (/-contour in A, and suppose that Z' m (A', z) has been defined by 
(4.9) for all m G S and all A' ^ A. Let us further assume by induction that Z' q (A', z) ^ for all 
m G S and all A' ^ A. We then define a smoothed cutoff function 4> q (Y, z) by 

4> q (Y,z)=l\ Xq ;m(Y,z), (4.10) 



where 



\< r .m{Y. :)-■ \ ( T - + T^ l °S 



Z'(IntY,z)e q (z)\Y\ 

z' m (intY,z)e m {z)\y\\) ■ v*- ll > 



Here Xq;m{Y, z) is interpreted as 1 if 9 m (z) or Z' m (lnt Y, z) is zero. 

As a consequence of the above definitions and the fact that Int m Y ^ A for all m G S, the 
expressions 

K>(Y,z) = Pz (Y)e q (z)-^4> q (y^) I! Z ^ m Y Z l (412) 

and 

iv~'(Y, 2;) = < ' ^ tt|A ^'^- e ' (4.13) 

10, otherwise, 

are meaningful for all z with q (z) / 0. By Lemma 3.9 we now know that Z' q (A, z) 7^ and the 
inductive definition can proceed. 

In the exceptional case 6 q (z) = 0, we let K' q (-, z) = K' q (-, z) = and z) = 0. Note that 
this is consistent with 4> q (Y, z) = 0. 

Remark 6. Theorem 4.2 stated and proved below will ensure that \K' q {Y,z)\ < e -(co+r/2)\Y\ for 
all ^-contours Y and all q <E S, provided r > Acq + 16. Hence, as it turns out a posteriori, the 
second alternative in (4.13) never occurs and, once we are done with the proof of Theorem 4.2, we 
can safely replace K q everywhere by K' q . The additional truncation allows us to define and use 
the relevant metastable free energies before stating and proving the (rather involved) Theorem 4.2. 
An alternative strategy would be to define scale dependent free energies as was done e.g. in [6]. 

4.2 Metastable free energies. 

Let us rewrite Z' q (A, z) as 

Z' q (A,z)=e q (z)^Z' q (A,z) (4.14) 

where 

z' q (A, z )= Yl WK^ z )- (4- 15 ) 

YeC(A,g) Y&Y 

We then define 

( q (z) = e q (z)e s ^ z \ (4.16) 
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where 

s q (z)= lim ±-AogZ' q (A,z) (4.17) 

|ah^,^oI a I 

By Lemma 3.10, the partition functions Z' q (A, z) and the polymer pressure s q (z) can be analyzed 
by a convergent cluster expansion, leading to the following lemma. 

Lemma 4.1 For each q G S and each z G G, the van Hove limit (4.17) exists and obeys the 
bound 

\s q (z)\ < e~ T / 2 . (4.18) 
If A is a finite subset ofL d and 9 q (z) / 0, we further have that Z' q (A, z) / and 

|log(C,(zr |A| ^(A,z))| < e-^ 2 |SA|, (4.19) 
while Cq(z) = and Z' q (A, z) = ifO q (z) = 0. 

Proof. Recalling the definition of compatibility between (/-contours from the paragraph before 
Lemma 3.9, C(A, q) is exactly the set of all compatible collections of (/-contours in A. Using 
the bound (4.13), the statements of the lemma are now direct consequences of Lemma 3.10, the 
definition (4.16), the representation (4.14) for Z' q (A, z) and the fact that we set K' q (Y, z) = if 
9 q (z) = 0. ' ' □ 

The logarithm of ( q {z) — or at least its real part — has a natural interpretation as the metastable 
free energy of the ground state q. To state our next theorem, we actually need to define these (and 
some other) quantities explicitly: For each z G G and each q G S with 9 q (z) / 0, let 

/,(*) = - log|C,W|, 

/(*) = min/ m (z), (4 _20) 
a q (z) = f q (z) - f(z). 

If 9 q (z) = 0, we set f q (z) = oo and a q = oo. (Note that sup zg ^ f(z) < oo by (4.16), the bound 
(4.18) and our assumption that 9{z) = max,, is bounded away from zero.) 

In accord with our previous definition, a phase q is stable at z if a q (z) = 0. We will also say 
that a (/-contour Y is stable at z if K' q (Y, z) = K q (Y, z). As we will see, stability of the phase q 
implies that all (/-contours are stable. Now we can formulate an analogue of Theorem 3.1 of [5] 
and Theorem 1.7 of [21]. 

Theorem 4.2 Suppose that r > 4co + 16 where cq is the constant from Lemma 3.9, and let 
e = e _T / 2 . Then the following holds for all z 6 G: 

(i) For all q G S and all q-contours Y, we have \K' q (Y, z)\ < e~( C0+T / 2 )l y l and, in particular, 

K' q (Y,z)=K' q (Y,z). 

(ii) IfY is a q-contour with a q (z) diam Y < j, then K' q (Y, z) = K q (Y, z). 

(iii) If a q (z) diam A < \, then Z q (A, z) = Z' q (A, z) ^ and 

\Z q {A,z)\ > e -/ w (*)|A|-g|0A|_ (421 ) 
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(iv) Ifm<ES, then 



\Z m {Kz)\ < e -/(-)|Al e 2^M. 



(4.22) 



Before proving Theorem 4.2, we state and prove the following simple lemma which will be 
used both in the proof of Theorem 4.2 and in the proof of Proposition 4.5 in the next subsection. 



Lemma 4.3 Let m, q G S, let z G & and let Y be a q-contour. 

(i) If(j) q (Y,z) > 0, then 

a 9 (|IntF| + |y|) < (r/4 + 2 + 4e- r/2 )|Y|. 

(ii) If(f> q (Y,z) > and Xq ;m (Y, z) < 1, then 

a m (|IntY| + |Y|) < (1 + 8 e - r/2 )|Y|. 



(4.23) 



(4.24) 



Proof of Lemma 4.3. By the definitions (4.10) and (4.11), the condition 4> q (Y, z) > implies that 

Z' n (lntY,z)e n {z)\ y \' 



max log 

nGS 



Z' q (lntY,z)e q (z)\ Y \ 



<(2 + r/4)m 



(4.25) 



Next we observe that (f) q (Y, z) > implies 9 q {z) / 0. Since the maximum in (4.25) is clearly 
attained for some n with 9 n (z) / 0, we may use the bound (4.19) to estimate the partition 
functions on the left hand side of (4.25). Combined with (4.16), (4.18), (4.20) and the estimate 
|<9Int Y| < \Y\, this immediately gives the bound (4.23). 
Next we use that the condition Xq;m{Y, z) < 1 implies that 



log 



Z!JkAY,z)9 ri 



\Y\ 



> (l + r/4)|Y| 



(4.26) 



Z g (lntY,z)e q (z)\ Y \ 

Since (4.26) is not consistent with 9 m (z) = 0, we may again use (4.19), (4.16), (4.18) and (4.20) 
to estimate the left hand side, leading to the bound 

(/ ? -/ m )(|IntY| + |Y|) > (r/4 + l-4e- T / 2 )|Y|. (4.27) 

Combining (4.27) with (4.23) and expressing a m as a q — (f q — f m ), one easily obtains the bound 
(4.24). ' ' □ 

As in [5], Theorem 4.2 is proved using induction on the diameter of A and Y. The initial 
step for the induction, namely, (i-ii) for diamY = 1 — which is trivially valid since no such 
contours exist — and (iii-iv) for diamA = 1, is established by an argument along the same lines 
as that which drives the induction, so we just need to prove the induction step. Let N > 1 and 
suppose that the claims (i-iv) have been established (or hold automatically) for all Y', A' with 
diam Y', diamA' < N. Throughout the proof we will omit the argument z in f m (z) and a m (z). 

The proof of the induction step is given in four parts: 
Proof of ( i). Let Y be such that diam Y = N. First we will show that the second alternative in 
(4.13) does not apply. By the bounds (4.1) and (4.18), we have that 



p z (Y)9 q (z)-\ Y \ 



< e 



-r\Y\ 



0(z) 
\0 q (z) 



\y\ 



< e -(r-2e)\Y\ e a q \Y\ 



(4.28) 
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while the inductive assumption (iv), the bound (4.19) and the fact that Y^ m I Int m Y\ = | Int Y\ 
and J2 m \dlnt m Y\ = \dhAY\ < \Y\, imply that 



n 



Z m (Int m Y,z) 



< e a q |Inty| e 36|Y|_ 



Assuming without loss of generality that <j> q (Y, z) > (otherwise there is nothing to prove), we 
now combine the bounds (4.28) and (4.29) with (4.23) and the fact that e = e~ r/2 < 2/r < 1/8, 
to conclude that \K' q (Y, z)\ < e -(t r -f- 5e ")l Y l < e -(| T - 4 )l y l. By the assumption r > 4c + 16, 
this is bounded by e~( C()+T / 2 )l y l, as desired. □ 

Proof of ( ii). Let diam Y = N and suppose that Y is a g-contour satisfying a q diam Y < r/4. 
Using the bounds (4.18) and (4.19), the definitions (4.16) and (4.20), and the fact that \dlnt Y\ < 
\Y\ we can conclude that 



max - — - log 

m£S \Y\ 



Z' m (lntY,z)e m (z) 



\y\ 



Z'(lntY,z)6 q (z)\ Y \ 



suppyuIntFl , r 



<»a ,„i +4t<^ + l. (4.30) 



In the last inequality, we used the bound | supp Y U Int Y\ < \Y\ diam Y, the assumption that 
a q di&mY < r/4 and the fact that 4e < 1. We also used that a q < oo implies 6 q / 0, which 
justifies the use of the bound (4.19). By the definitions (4.10) and (4.11), the bound (4.30) implies 
that <f) q (Int Y, z) = 1. On the other hand, Z q (lnt m Y, z) = Z' q (lnt m Y, z) for all m £ S by the 
inductive assumption (iii) and the fact that diam Int m Y < diam Y = N. Combined with the 
inductive assumption (i), we infer that K' q (Y, z) = K' q (Y, z) = K q (Y, z). □ 

Proof of (iii). Let A C 7h d be such that diam A = and a q diam A < r/4. By the fact that (ii) is 
known to hold for all contours Y with diam Y < N, we have that K' q (Y, z) = K q (Y, z) for all Y 
in A, implying that Z q (A, z) = Z' q (A, z). Invoking (4.19) and (4.20), the bound (4.21) follows 
directly. □ 

Proof of (iv). Let A be a subset of 7L d with diam A = N. Following [21, 5], we will apply the 
cluster expansion only to contours that are sufficiently suppressed and handle the other contours 
by a crude upper bound. Given a compatible collection of contours Y, recall that internal contours 
are those contained in Int Y of some other Y € Y while the others are external. Let us call 
an m-contour Y small if a m diam Y < r/4; otherwise we will call it large. The reason for this 
distinction is that if Y is small then it is automatically stable. 

Bearing in mind the above definitions, let us partition any collection of contours ¥ £ M(A, m) 
into three sets Y^UY^UYj^ of internal, small-external and large-external contours, respec- 
tively. Fixing Y^g e and resumming the remaining two families of contours, the partition function 
Z m (A, z) can be recast in the form 

Z m (A,z) = ^Zr all (Ext,z) I] [ Pz (Y] H Z n (Int n y)}. (4.31) 

Y YeY ™ e5 

Here the sum runs over all sets Y of mutually external large m-contours in A, the symbol Ext = 
ExtA(Y) denotes the set f|y e Y( Ext ^nA) and Z™ all (Ext, z) is the partition sum in Ext induced 
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by Y. Explicitly, Z^ a (A, z) is the quantity from (3.6) with the sum restricted to the collections 
Y G M(A, m) for which all external contours are small according to the above definition. 

In the special case where 9 m {z) = 0, all contours are large by definition (recall that a m = oo 
if 9 m vanishes) and the partition function Z™ all (A, z) is defined to be zero unless A = 0, in 
which case we set it to one. We will not pay special attention to the case 6 m = in the sequel of 
this proof, but as the reader may easily verify, all our estimates remain true in this case, and can 
be formally derived by considering the limit 

Using the inductive assumption (iv) to estimate the partition functions Z n (Int n Y), the Peierls 
condition (4.1) to bound the activities p z {Y), and the bound (4.18) to estimate 6{z) by e~^e e , 
we get 

II {p&) II Z n(^nY)} < n {e^le^^M^ 3 ^} 

YeY ne5 YeY 



-/|A\Ext| JT e -(r-3e)|Y| 



(4.32) 

YeY 

Next we will estimate the partition function Z™ all (Ext, z). Since all small m-contours are stable 
by the inductive hypothesis, this partition function can be analyzed by a convergent cluster ex- 
pansion. Let us consider the ratio of Z™ all (Ext, z) and Z^(Ext, z). Expressing the logarithm of 
this ratio as a sum over clusters we obtain a sum over clusters that contain at least one contour of 
size \Y\ > diam Y > r/a m > 2/a m . Using the bound (3.21) with r\ = r/2 we conclude that 

Z™ all (Ext,z) 



< e |Ext|e-/""\ (4.33) 



Z^Ext,*) 

Combined with Lemma 4.1 and the definitions (4.20), this gives 

|^ mall (Ext,z)| < e -(/ m - e - T/a -)|Ext| e£ ~|9A| JT e e|Y|_ (4 _ 34) 

YeY 

We thus conclude that the left hand side of (4.31) is bounded by 



\Z m (A,z)\ <max(V^/ 2 )l Ext l [] e^'^ Y \\ 
Y YeY 

^ e -6|Ext| Yl 



x e -/|A| e e|SA| ^ e -6|Ext| T"f ( -(or/.l l< )\) 

Y YeY 



(4.35) 



where b = a m /2 - e T l am . Note that b > e T l am which is implied by the fact that 4e T l am < 
Aa m /r < a m . 

For the purposes of this proof, it suffices to bound the first factor in (4.35) by 1. In a later 
proof, however, we will use a more subtle bound. To bound the second factor, we will invoke 
Zahradnfk's method (see [21, Main Lemma] or [5, Lemma 3.2]): Consider the contour model with 
weights K(Y) = e -( 3r / 4 - 4 )l y l if Y is a large m-contour and K(Y) = otherwise. Let Z(A) be 
the corresponding polymer partition function in A — see (3.11) — and let <p be the corresponding 
free energy. Clearly Z(A) > 1 so that — ip > 0. Since 3r/4 — 4 > cq + r/2, we can use 
Lemmas 3.9 and 3.10 to obtain further bounds. For the free energy, this gives < —if < 



34 



M. BISKUP, C. BORGS, J.T. CHAYES, AND R. KOTECKY 



min{e, e T / a ™} because the weights of contours smaller than 2/a m identically vanish. Since 
b > e~ T l am , this allows us to bound the sum on the right hand side of (4.35) by 

^ eV |Ext| Yl e -(3r/4-4e)|y| < Ext | "Q |Y| e -(3r/4-5e)|Y| j (4 36) 

Y YGY Y YgY 

Using Lemma 3.10 once more, we have that Z(Inty)e ¥ ' |Intl V~ |y| > 1. Inserting into (4.36), 
we obtain 



^ e -6|Ext| -Q e -(3r/4-4e)|Y| < J- ^(| Ext |+E ye5 (I Int Y| + | Y|)) JJ j ^ Y )K{Y) } 

ygy y ygy 



Y YgY 

(4.37) 

Consider, on the other hand, the polymer partition function Z{A) in the representation (3.11). 
Resuming all contours but the external ones, we obtain precisely the right hand side of (4.37), 
except for the factor e^l A L This shows that the right hands side of (4.37) is equal to Z(A)e (p \ A \ 
which — again by Lemma 3.10 — is bounded by e e l 9A l Putting this and (4.35) together we obtain 
the proof of the claim (iv). □ 

4.3 Differentiability of free energies. 

Our next item of concern will be the existence of two continuous and bounded derivatives of the 
metastable free energies. To this end, we first prove the following proposition, which establishes 
a bound of the form (4.22) for the derivatives of the partition functions Z m (A, z). 

Proposition 4.4 Let r and M be the constants from (4.1) and (4.3), let e = e~ T / 2 , and suppose 
that t > 4co + 16 where cq is the constant from Lemma 3.9. Then 

didlz m {A,z)\ < e-^)l A l(2M|A|) W >1 aA l, (4.38) 
holds for all z G 6, all m G S, and all £, t > with I + t < 2. 

Proof. Again, we proceed by induction on the diameter of A. We start from the representation 
(4.4) which we rewrite as 

z m {Kz)=Y, n m*o n z ( y > 2 ^ < 4 - 39 ) 

Y ext zGExt YGY exl 

where we abbreviated Z(Y, z) = p z (Y) Y[ n Z n (lnt n Y, z). Let 1 < t < oo be fixed (later, we 
will use that actually, i < 2) and let us consider the impact of applying d\ on Z m (A, z). Clearly, 
each of the derivatives acts either on some of m 's, or on some of the Z(Y, z)'s. Let k x be the 
number of times the term 9 m (z) is differentiated "at x" and let iy be the number of times the 
factor Z(Y,z) is differentiated. Let k = (k x ) and i = (iy) be the corresponding multiindices. 
The resummation of all contours Y for which iy = and k x = for all x 6 supp Y U Int Y 
then contributes a factor Z m (Ext\(Y ex )\A',z), where we used Y ex to denote the set of all those 
Y G Y ext for which iy > 0, Ext A (Y ext ) = A\ \J Yer » (supp Y Ulnt Y), and A' = {x: k x > 0}. 
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(Remember the requirement that no contour in ExtA(Y ) \ A' surrounds any of the "holes.") 
Using this notation, the result of differentiating can be concisely written as 

d e z Z m (^,z) = J2 E Z m (Ext A (Y ext )\A',z) 

Y 8 *' A'CExt A (Y cxl ) 

X E fcT-T II^-M*) II 9rZ(Y,z). (4.40) 

k,i x£A' YG¥ ex ' 

k+i=i 

Here the first sum goes over all collections (including the empty one) Y 6 "' of mutually external 
contours in A and the third sum goes over all pairs of multiindices (k, i), k x = 1, 2, . . . , x G A', 
iy = 1,2,..., Y G Y CXt . (The terms with |A'| + IY 6 "'] > ^ vanish.) We write k + i = I 
to abbreviate Yl x k x + ^ y iy = I and use the symbols fc! and i\ to denote the multi-index 
factorials \\ x kx* and Y\ Y ^yU respectively. 

We now use (4.3) and (4.18) to bound \d%°O m (z)\ by (M) kx e l e~^ z \ Employing (4.1) and 
(4.18) to bound the derivatives of p z (Y), and the inductive hypothesis to bound the derivatives of 
Z m (Int m Y,z), we estimate \d^Z(Y,z)\ by [2M|y(y)|]^ e -( r - 3e ")l Y l e -^ 2 )l y ( Y )l (recall that 
V(Y) was defined as supp Y U Int Y). Finally, we may use the bound (4.22) to estimate 

|i? m (Ext A (r xt )\A^)|<e 2 ^^ (4.41) 
Combining these estimates and invoking the inequality 

|5(Ext A (r xt )\A , )|<|5A| + |A / |+ Y, l y l> ( 4 - 42 ) 

ye¥ cxt 

we get 

^ ro (A,,)|< e *1»l e -/MI*l£ y, 

Y 6 "' A'cExt A (Y eXt ) fe fe ^ 

x ]J(Me 3 ~ e ) k * 11 (2M\V{Y)\) iY e- {T -^ Y \. (4.43) 

xeA' YG¥ ex ' 

Let us now consider the case i = 1 and 1 = 2. For t = 1, the sum on the right hand side of 
(4.43) can be rewritten as 

^(Me 3? + 2Me-( r - 5? )l y l), (4.44) 

xeA Y:xeV(Y)cA 

while for i = 2, it becomes 

Y ({Me 3i ) 2 + 2Me 3i 2M £ e -0-K)IH + ( 2M ) 2 £ J] e -(^")l y l) , (4.45) 

yeV(Y)cA 

where the last sum goes over sets of mutually external contours Y 6 " in A such that {x, y} C 
U ye ^ext V(Y) and {x, y} n V(Y) / for each Y G Y ext . Note that the last condition can only 
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— ext 

be satisfied if Y contains either one or two contours. Introducing the shorthand 

S= e-^- 5? )l y l (4.46) 

Y:0eV(Y)cZ d 

we bound the expression (4.44) by (e u + 2S')M|A|, and the expression (4.45) by (e 6i + Ae u S + 
4(S + S 2 ))M 2 \A\ 2 . Recalling that c was defined in such a way that the bound (3.28) holds, we 
may now use the fact that r — 5e — cq > \t to bound S by e _2 e. Since e < 1/8, this implies that 
the above two terms can be estimated by (e 3/8 + \e~ 2 )M\A\ < 2M\A\ and (e 6/8 + ie 3/8 ~ 2 + 
\{e- 2 + ie~ 4 ))M 2 |A| 2 < 4M 2 |A| 2 , as desired. 

This completes the proof for the derivatives with respect to z. The proof for the derivatives 
with respect to z and the mixed derivatives is completely analogous and is left to the reader. □ 

Next we will establish a bound on the first two derivatives of the contour weights K'. Be- 
fore formulating the next proposition, we recall the definitions of the polymer partition function 
Z' q (A, z) and the polymer pressure s q in (4.17) and (4.15) . 

Proposition 4.5 Let r and M be the constants from (4.1) and (4.3), let Co be the constant from 
Lemma 3.9, and let e = e~ T ^ 2 . Then there exists a finite constant T\ > Acq + 16 depending 
only on M, d and \S\ such that if t > t\, the contour weights K' q (Y, ■) are twice continuously 

differentiate in 0. Furthermore, the bounds 

\d l z dlK' q (Y,z)\ < e -(-o+-/2)i^i (4 _ 47) 

and 

\di4K(^,z)\ < |A|^>WI A l+^ A l (4.48) 
hold for all q G S, all z G G, all q-contours Y, all finite AcZ 1 * and all £, i > with £ + 1 < 2. 

Proposition 4.5 immediately implies that the polymer pressures s q are twice continuously dif- 
ferentiable and obey the bounds of Lemma 3.10. For future reference, we state this in the follow- 
ing corollary. 

Corollary 4.6 Let t\ be as in Proposition 4.5. If 't > t\ and q G S, then s q is a twice continu- 
ously differentiable function in & and obeys the bounds 

\d w s q \<e- r/2 and \d w d w ,s q \ < e~ r/2 , (4.49) 

valid for any w, w' G {z, z} and any z G G. 

Proof of Proposition 4.5. Let r > t\ > 4co + 16. Then Theorem 4.2 is at our disposal. It will be 
convenient to cover the set & by the open sets 

= {z^G: \Q q {z)\ < e-( r / 4+2+6e "^(z)} (4.50) 

and 

0^ = {z£S': \9 q (z)\ > e-( r / 4+2+8e "^(z)}. (4.51) 

We first note that K' q (Y, z) = if z G . Indeed, assuming K' q (Y, z) / we necessarily have 
4> q {Y, z) > 0, which, by (4.23), implies that a q < r/4 + 2 + 4e and thus log 9(z) - log \9 q (z)\ < 
r/4 + 2 + 6e, which is incompatible with z G 6^ . Hence, the claims trivially hold in and 
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it remains to prove that K' q (Y, •) is twice continuously differentiable in and that (4.47) and 

(4.48) hold for all z G d^f 1 . As in the proof of Theorem 4.2 we will proceed by induction on the 

diameter of Y and A. Let N > 1 and suppose that K' q (Y, •) G C 2 (^ 9) ) and obeys the bounds 
(4.47) for all q G S and all g-contours Y with diam Y < N, and that (4.48) holds for all q G 5 
and all A C Z d with diam A < N - 1. 

We start by proving that K' (Y, •) G C 2 {^^) whenever Y is a (/-contour Y of diameter N. To 

this end, we first observe that in we have that 9 q (z) / and hence also Z' q (lnt Y, z) / 0. 
Using the inductive assumption, this implies that the quotient 

Z' m (lntY,z)e m (z)\ Y \ 
Qm > Y{2) ~ Z>QntY t z)O q (z)\y\ (4 " 52) 

is twice continuously differentiable in S 1 ^, which in turn implies that Xq;m{Y, z) is twice con- 
tinuously differentiable. Combined with the corresponding continuous differentiability of p z (Y), 
9 q {z), Z m (lnt m Y, z), and Z' (Int m Y, z), this proves the existence of two continuous derivatives 
of z i — ^ K' q (Y, z) with respect to both z and z. 

Next we prove the bound (4.48) for diam A = iV — 1. As we will see, these bounds follow im- 
mediately from the inductive assumptions (4.47) and Lemma 3.10. Indeed, let lq{Y) = K' q (Y, z) 
if diamY < N — 1, and l q (Y) = if diamY > N — 1. The inductive assumptions (4.47) 
then guarantee the conditions (3.22) of Lemma 3.10. Combining the representation (3.15) for 
log^(A, z) with the estimate (3.23) from Lemma 3.10 we thus conclude that 

\didl\ogZ' q (k,z)\ < |A|e, (4.53) 

while (3.26) gives the bound 

\Z' q {A,z)\ <e s " W+il9M . (4.54) 

Combining these bounds with the estimates e|A| < |A| and e 2 1 A | 2 + e|A| < |A| 2 , we obtain the 
desired bounds (4.48). 

Before turning to the proof of (4.47) we will show that for z G ff^, the bound (4.48) implies 

\d f z dlz>(A,z)\ < (M ie -/ 4+3 |A|) £+£ e-^( 2 )l A l +? l 9A l (4.55) 

with M\ = 1 + M. Indeed, invoking the assumption (4.3), the definition of &^ , and the fact 
that e < 1/8, we may estimate the first and second derivative of ^(z)l A l by 

didle q (zr\\ < (mM^y^ IWI |A| < (M\A\e^) W \9 q (z)\M. (4.56) 

Combined with (4.14) and (4.48) this gives (4.55). 

Let Y be a g-contour with diamY = N, and let us consider the derivatives with respect to z; 
the other derivatives are handled analogously. By the assumption (4.1) and the bound (4.18), we 
have 

\dl Pz (Y)\ < |Y|^e-^- 2? )l y le^l y l|^(z)|l y l, (4.57) 
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while (4.3) and the assumption that z G fi^ ( cf ( 4 - 56 )) yields 

\die q {z)-\ Y \ \ < {\Y\ + l) e (Me T ^ +3 )%(z)\~ lYl . (4.58) 
Further, combining the bound (4.55) with Theorem 4.2 and Proposition 4.4 we have 

% TT Z ^ lntmY ' z) < |Inty| £ (2M + 2M ie 2? l y le 3 +^ 4 )V ? l y le^l Inty l. (4.59) 
meS Z q {lnt m Y,z) 

Finally, let us consider one of the factors x g;m (Y, z). To bound its derivative, we may assume 
that z is an accumulation point of z' with Xq;m{Y, z') < 1 (otherwise its derivative is zero), so by 
Lemma 4.3(h) we have that a m < 1 + 8e and thus log 6{z)- log \9 m (z)\ < l + 10e < r/4+2+8e, 
implying that z G . We may therefore use the bounds (4.56) and (4.55) to estimate the 
derivatives of Xq-,m(X, z), yielding the bound 

\dix q -, m {Y,z)\ < C{\ IntY| + \Y\Y (4M ie 3+r / 4 e 2?|y| )^ (4.60) 



where C is a constant bounding both the first and the second derivative of the mollifier function x- 
Combining all these estimates, we obtain a bound of the form 

\d e z K' q (Y,z)\ < C(\ IntY| + |y|)^W4 e -(r-c e -)|y| e a 9 (|IntY|+|Y|) {4M) 

with a constant C that depends on M and the number of spin states |<S|, and a constant c that 
depends only on |<S|. Using the bound (4.23) and the fact that e £r/4 < e (r/8)|y| (note that \Y\ > 
(2R + l) d > 4 by our definition of contours), we conclude that 

\diK' q (Y,z)\ < C(|Inty| + |y|)« e -(5r/8-3-c6)|y|_ (4 _ 62) 

Increasing t\ if necessary to absorb all of the prefactors, the bound (4.47) follows. □ 

We close the subsection with a lemma concerning the Lipschitz continuity of real- valued func- 
tions zh j[z) and z i-> e^ aq ^ on &: 

Lemma 4.7 Let t\ be as in Proposition 4.5 and let M\ = AM + 1. If r > r±, q G S, and if 
z,zo G & are such that [zq, z] = {sz + (1 — s)z$ : < s < 1} C G, then 

\f(z )-f(z)\<M 1 \z-z \ (4.63) 

and 

| e -««W - e -««(«o)| < 2M x \z - zo | e Ml|z_Zo1 . (4.64) 

Proof Let Cq( z ) De tne quantity defined in (4.16), and let e = e~ T l 2 . Combining the assumption 
(4.3) with the bounds (4.49) and (4.18), we get the estimate 

\d w ( q {z)\ < {Me 2i + e)e- f ( z \ w, w' G {z, z}. (4.65) 
With the help of the bound Me 2i + e < 2M + 1/2 = Mi/2, we conclude that 



| e -/«(*0 _ e -/ 9 (^)| < Mi / e" /(z,) Id/I, «i,«2 G [z ,z], (4.66) 

■'[21,22] 
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where \dz'\ denotes the Lebesgue measure on the interval [z , z\. Using that / = max, f q , this 
implies 

| e -/(*i) _ e - f( - Z2) \ <M X \ e~ f{z ' ] \dz'\, z u z 2 G [zo, z]. (4.67) 

Now if (4.63) is violated, i.e., when \f(z) — f(zo)\ > (Mi + e)\z — zq\, then the same is true 
either about the first or the second half of the segment [zq , z] . This shows that there is a sequence 

of intervals [zi, n ,z 2 , n ] of length 2~ n |z - z\ where \f(z hn ) - f(z 2 , n )\ > {Mi + e)\z 1:Tl - z 2 , n \- 
But that would be in contradiction with (4.67) which implies that 

lim i/(;,->-/( r >i _ lim < Mi, ,4^ 8 ) 

n-oo \z ltn -Z2, n \ ™ J[* W2 , n ] e ~ /( Wl 

where we use the mean-value Theorem and a compactness argument to infer the first equality. 
Hence, (4.63) must be true after all. 

To prove (4.64), we combine the triangle inequality and the bound f q (zo) > f(zo) with (4.66) 
and (4.67) to conclude that 

| e -a,(z) _ e -a,(z )| _ |g/W e -/q(^) _ e /(zo) e -/ 9 (zo)| 

< e /W| e -/,W _ e -/,(*o)i + l[ q{Z °l . \ e -fM - e~f^\ r4 m 

1 1 e -f(z) e -f( z o) 1 1 K^.W) 

< 2Mi e f ^- f{z ">\dz'\. 

Jz 

Bounding f(z) — f(z') by M\\z — zq\, we obtain the bound (4.64). □ 
4.4 Torus partition functions. 

In this subsection we consider the partition functions Z q (A,z), defined for A C Tl in (3.6). 
Since all contours contributing to Z q (A,z) have diameter strictly less than L/2, the partition 
function Z q (A, z) can be represented in the form (4.8), with K q (Y, z) defined by embedding the 
contour Y into Z d . Let Z' (A, z) be the corresponding truncated partition function, defined with 
weights K' q (Y, z) given by (4. 12). Notice, however, that even though every contour Y C A can be 
individually embedded into the relation of incompatibility is formulated on torus. The poly- 
mer partition function Z' q {A, z) and Z' (A, z) can then again be analyzed by a convergent cluster 
expansion, bearing in mind, however, the torus incompatibility relation. The torus analogue of 
Lemma 4.1 is then as follows: 

Lemma 4.8 Assume that r > t\, where t\ is the constant from Proposition 4.5 and let q G S 
and z G 6 he such that 9 q (z) / 0. Then 

4 log (( q (zy^Z' q (A,z)) I < e^/ 2 |<9A| + 2|A|e-^ 4 (4.70) 
for any A C T^, any z £ 6, t = 0, 1, and w G {z, z}. 

Proof. Let us write Z'JA,z) in the form (4.14). Taking into account the torus compatibility 
relation when comparing the cluster expansion for logZ' q (A, z) with the corresponding terms 
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contributing to s q \ A|, we see that the difference stems not only from clusters passing through the 
boundary dA, but also from the clusters that are wrapped around the torus in the former as well 
as the clusters that cannot be placed on the torus in the latter. For such clusters, however, we 
necessarily have J2y XOOI^I — Since the functional $(Y) = K' q (Y, z) satisfies the bound 
(3.20) with 77 = t/2, we may use the bound (3.21) to estimate the contribution of these clusters. 
This yields 

log^(A,z)-a,|A|| < e - r / 2 |9A|+2|A|e- rL / 4 , (4.71) 

which is (4.70) for t = 0. To handle the case I = 1, we just need to recall that, by Proposition 4.5, 
the functional }(Y) = K' q {Y,z) satisfies the bounds (3.22) with 77 = t/2. Then the desired 
estimate for i = 1 follows with help of (3.23) by a straightforward generalization of the above 
proof of (4.71). □ 

Next we provide the corresponding extension of Theorem 4.2 to the torus: 

Theorem 4.9 Let r > 4co + 16 where c$ is the constant from Lemma 3.9, and let us abbreviate 
e = e~ T l 2 . For all z G G, the following holds for all subsets A of the torus T^: 

(i) Ifa q (z) diam A < |, then Z q (A, z) = Z' q (A, z) ^ and 



(ii) IfmeS, then 



\Z q (A,z)\ > e -/,W|A| e - e -|aA|-2|A|e-V4_ (4 _ 72) 

\Z m (A, Z)\ < e -/W|A|+2 e -|aA|+4|A|e-V4_ {4 J?>) 

(iii) IfmeS, then 

\Z m (T L , z)\ < e~^ Ld max{ e -^W Ld / 2 , e -^/AyL^-^_ {4J4) 

Remark 7. The bounds (4.72) and (4.73) are obvious generalizations of the corresponding bounds 
in Theorem 4.2 to the torus. But unlike in Proposition 4.5, we will not need to prove the bounds 
for the derivatives with respect to z. When such bounds will be needed in the next section, we 
will invoke analyticity in z and estimate the derivatives using Cauchy's Theorem. 

Proof of (i). Since all contours can by definition be embedded into Z d , Theorem 4.2(h) guarantees 
that K' q (Y, z) = K q (Y, z) for all g-contours in A and hence Z q (A, z) = Z' q (A, z). Then (4.72) 
follows by Lemma 4.8 and the definition of f q . □ 

Proof of (ii). We will only indicate the changes relative to the proof of part (iv) of Theorem 4.2. 
First, since all contours can be embedded into 7L d , we have that a corresponding bound — namely, 
(4.22) — holds for the interiors of all contours in A. This means that all of the derivation (4.31— 
4.35) carries over, with the exception of the factor e e ' aA ' in (4.34) and (4.35) which by Lemma 4.8 
should now be replaced by e e l i9A l+ 2 l A l e tL/4 . i n order to estimate the last sum in (4.35), we will 
again invoke the trick described in (4.36-4.37). This brings in yet another factor e e l 9A l+ 2 l A l e tL/a _ 
From here (4.73) follows. □ 

Proof of (iii). The estimate is analogous to that in (ii); the only difference is that now we have to 
make use of the extra decay from the maximum in (4.35). (Note that for A = we have |<9A| = 
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and I A| = L d .) Following [5], this is done as follows: If Y is a contour, a standard isoperimetric 
inequality yields 

1*1 > ^i\0(snppY UlntY)| > | suppF U Int Y\^r . (4.75) 

Hence, if Y is a collection of external contours in Tl and Ext is the corresponding exterior set, 
we have 

d-l 

\ Y \ ^ |suppyulnty|^ > (j2 |suppyulnty| J d = (L d -| Ext l)^- (4.76) 

YeY YGY ygy 

Writing | Ext | = (1 — x)L d where x 6 [0, 1], the maximum in (4.35) is bounded by 

sup exp {-^L d (l _ x) _ T-L^x^r). (4.77) 

X6[0,l] 12 4 J 

The function in the exponent is convex and the supremum is thus clearly dominated by the bigger 
of the values at x = and x = 1. This gives the maximum in (4.74). □ 

Apart from the partition functions Z m (TL, z), we will also need to deal with the situations 
where there is a non-trivial contour network. To this end, we need a suitable estimate on the 
difference 

Z h ^(z) = Zf\z) - Z m (T L , z). (4.78) 
This is the content of the last lemma of this section. 

Lemma 4.10 There exists a constant cq depending only on d and \S\ such that for r > 4co + 16 

and all z G 6, we have 

\Z h l\z)\ < L d e-^e 5Lde ' TL/4 C(z) Ld . (4.79) 
Proof. Let cq be the constant from Lemma 3.9, and let cq = co(d, \S\) > Co be such that 

(|5|e" C0 )l A l < L d , (4.80) 

AcT L 

where the sum goes over all connected subsets A of the torus T l (the existence of such a constant 
follows immediately from the fact that the number of connected subsets A C Z d that contain a 
given point x and have size k is bounded by a <i-dependent constant raised to the power k). 

The proof of the lemma is now a straightforward corollary of Theorem 4.9. Indeed, invoking 
the representation (3.8) we have 

Z h l g {z)= Pz(N) Hz m (A m ((&,N),z), (4.81) 

(%,JS)eM L rneS 
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where A m (0, N) is defined before Proposition 3.7. Using (4.2) and (4.73) in conjunction with the 
bounds 9(z) < ((z)e 2i and £ me5 |<9A m (0, N)| < we get 

I^WI < CW Ld e 4Lde " Ti/4 ^ e -(T-4£)|>f|_ (4 _ 82) 

(«,N)eMi 

Taking into account that each connected component of supp ]Sf has size at least L/2, the last sum 
can be bounded by 

£ e -(r-4g)|>f| < ^ i-S" < 5e 5 (4.83) 
(0,W)GA4 L n=l ' 

where 

5= ^ (|S|e-(-- 4? )) |A| (4.84) 

AcT L 
|A|>L/2 

is a sum over connected sets A C Tl of size at least L/2. Extracting a factor e~ ri//4 from the 
right hand side of (4.84), observing that r/2 — 4e > Co, and recalling that cq was defined in such 
a way that (4.80) holds, we get the estimate S < L d e" rL/A . Combined with (4.82) and (4.83) 
this gives the desired bound (4.79). □ 



5. Proofs of main results 

We are finally in a position to prove our main results. Unlike in Section 4, all of the deriva- 
tions will assume the validity of Assumption C. Note that the assumptions (4.1-4.3) follow from 
Assumptions C0-C2, so all results from Section 4 are at our disposal. Note also that p z (Y), 
p z ("N) and 9 m {z) are analytic functions of z by Lemma 3.6, implying that the partition functions 
Z m (A, •) and Z^ CI are analytic functions of z. 
We will prove Theorems A and B for 

r = max{ri,4c + 16, 21og(2/a)} (5.1) 

where t\ is the constant from Proposition 4.5, Co is the constant from Lemma 4.10 and a is the 
constant from Assumption C. Recall that n > 4co + 16, so for r > To we can use all results of 
Section 4. 

First, we will attend to the proof of Theorem A: 

Proof of Theorem A. Most of the required properties have already been established. Indeed, let ( q 
be as defined in (4.16). Then (2.9) is exactly (4.18) which proves part (1) of the Theorem A. 

In order to prove that dz( q (z) = whenever z G ,5^ q , we recall that ( q (z) = 6 q {z)e Sq ^ 
where 9 q (z) is holomorphic in & and s q (z) is given in terms of its Taylor expansion in the contour 
activities K' (Y, z). Now, if a q (z) = — which is implied by z G 5^ q — then K' (Y, z) = K q (Y, z) 
for any ^-contour Y by Theorem 4.2. But dzK q (Y, z) = by the fact that p z (Y), Z q (lnt m Y, z) 
and Z m (lnt m Y, z) are holomorphic and Z q (lnt m Y, z) / 0. Since s q is given in terms of an 
absolutely converging power series in the K q 's, we thus also have that d z e Sq ( z } = 0. Hence 
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To prove part (3), let z G 5P m fl 5? n for some distinct indices m, n G TZ. Using Lemma 4.1 we 
then have 

e m (z) > e(z)e- 2e ~ T/2 (5.2) 

and similarly for n. Since a > 2e _r °/ 2 > 2e~ T / 2 , we thus have z G J§? a (m) n Jz? a (n). Using the 
first bound in (4.49), we further have 



Cm{z) Cn(z) 



> \d z e m {z) - d z e n {z)\ - 2 e - T / 2 . (5.3) 



Applying Assumption C3, the right hand side is not less than a — 2e T l 2 . Part (4) is proved 
analogously; we leave the details to the reader. □ 

Before proving Theorem B, we prove the following lemma. 
Lemma 5.1 Let e > 0, let t\ be the constant from Proposition 4.5, and let 

sf\z) = -^\ogZ' q {Y L ,z) (5.4) 



and 



C[ L \z) = e q {z)e^ L) ^. (5.5) 
Then there exists a constant M depending only on e and M such that 

\d"4 L) {z)\ < (e\) 2 (M y\c( L \z)\ (5.6) 

holds for all q G S, all £ > 1, all r > T\, all L > r /2 and all z G G with a q {z) < r/ (4L) and 
dist(z, > e. 

Proof. We will prove the lemma withe the help of Cauchy's theorem. Starting with the derivatives 
of 6 q , let e = min{e, 1 / (4Mi)} where = 1 + 4M is the constant from Lemma 4.7, and let 
2' be a point in the disc B> eo (z) of radius eo around z. Using the bounds (4.18) and (4.63), we 
now bound 

\0q{z')\ < e^ /(2,) < e t+Mico e -f(z) < e e+M 1 e +a q (z) e ^f q (z) < |^(^)| e 2e+S 1 e +a,(z)_ (5 _ 7) 

With the help of Cauchy's theorem and the estimates e < 1/8, Mie < 1/4 and a ? (z) < 1/2, this 
implies 

\d%(z)\ 



< lltfe 1 ' 4 * 1 ' 4 * 1 ' 2 < ei&e^y. (5.8) 



In order to bound the derivatives of , let us consider a multiindex X contributing to the cluster 
expansion of Sq L \ and let k = maxy. X (y) >0 diam Y. Defining 

e k = min{e, (20eMiA;)~ 1 }, (5.9) 

where M\ = 1 + AM is the constant from Lemma 4.7, we will show that the weight K' q (Y, •) of 
any contour Y with X(Y) > is analytic inside the disc D efc (z) of radius about z. Indeed, let 
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\z — z'\ < Combining the assumption a q (z) < r/(4L) < 1/2 with Lemma 4.7, we have 

e -a q {z>) > e -o,(«) _ 2eMie k > 1 - a q (z) - 2eM iek 

> 1 - ^max{a g (z),10eM iefc } > e -2max{a 9 ( 2 ),10eM iefc }_ 


Here we used the fact that x + y < | max{x, 5y} whenever x, y > in the last but one step, and 
the fact that e~ 2x < 1 — (1 — e _1 )2x < 1 — |x whenever x < 1/2 in the last step. We thus have 
proven that 

a q (z') < max{2a q {z),20eMie k } < max{^,i} < (5.11) 

so by Theorem 4.2, K' q (Y,z') = K q (Y,z') and Z q (Int m Y, z') / for all m G S and z' G 
D efe (z). As a consequence, K' q (Y, •) is analytic inside the disc D efe (z), as claimed. 

At this point, the proof of the lemma is an easy exercise. Indeed, combining Cauchy's theorem 
with the bound \K' q (Y,z')\ < e -( T / 2 + c «)l y l < e - c °l Y le-( T / 2 ) diamY , we get the estimate 

3iJ{K ! q (Y,z')^ < eieil[e-^ + ^\ Y W Y ) < ^e^^ k \{e~^ Y \ (5.12) 

Y Y Y 

Bounding e~ l eT^ l 2 ^ k by e^[ e k e e' k < (ie^e^Y, we conclude that 

\di \{K' q {Y,z' f {Y) \ < £!(£e- 1 er 1 )'n e_C ° |y|X(y) - ^ 13 ) 

Y Y 

Inserted into the cluster expansion for s q L ^ , this gives the bound 



disfXz) <i\{ie-\^)\ (5.14) 



which in turn implies that 



Combining this bound with the bound (5.8), we obtain the bound (5.6) with a constant Mq that 
depends only on e and M\ , and hence only on e and M. □ 

Next we will prove Theorem B. Recall the definitions of the sets y e (m) and ^4(Q) from 
(2.13) and (2.14) and the fact that in Theorem B, we set k = r/4. 

Proof of Theorem B(l-3). Part (1) is a trivial consequence of the fact that d m (z), /> 2 (N) and p z {Y) 
are analytic functions of z throughout &. 

In order to prove part (2), we note that z G ((/) implies that a q {z) < k/L = r/(4L) 
and hence by Theorem 4.2(h) we have that K' q (Y, z) = K q (Y, z) for any (/-contour contributing 

to Z q (TL, z). This immediately implies that the functions s q L ^ and ( q L \z) defined in (5.4) and 
(5.5) are analytic function in S fi K i L {q). Next we observe that r > 4co + 16 implies that tL/8 > 
t/8 > log 4 and hence 4e~ rL / 4 < e~ TL / 8 . Since z G ^ k /l( ( ?) implies a q (z) < oo and hence 
9 q {z) / 0, the bounds (2.15-2.16) are then direct consequences of Lemma 4.8 and the fact 
thatOT L = 0. 

The bound (2.17) in part (3) finally is nothing but the bound (5.6) from Lemma 5.1, while he 
bound (2.18) is proved exactly as for Theorem A. Note that so far, we only have used that r > tq, 
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except for the proof of (2.17), which through the conditions from Lemma 5.1 requires L > r/2, 
and give a constant Mq depending on e and M. □ 

Proof of Theorem B(4). We will again rely on analyticity and Cauchy's Theorem. Let Q C 1Z 
and let Q' C S be the set of corresponding interchangeable spin states. Clearly, if m and n are 

interchangeable, then (ffl = C,n^ and, recalling that q m denotes the set of spins corresponding 
to m G 1Z, we have 

S Q {z) = 2%*{z) ~ £ [Ck L) (z)] Ld = Zf\z) - £ Z' a (T L , z). (5.16) 

neQ' n£Q' 

Pick a zq G ^ k /l{Q)- F° r n £ Q'» we then have a n (zo) < r/ (4L), and by the argument leading 
to (5.11) we have that a n (z) < r/(2L) provided r/(4L) < 1/2 and 2eM x \z - z \ < \^ L . 
On the other hand, if m G S \ Q', then a m (zo) > rj (8L), and by a similar argument, we get 
that a m (z) > r/(16L) if r/(8L) < 1 and 2eM 1 \z - z \ < jq-^. Noting that r > f implies 
r > 4co + 16 > 16, we now set 

e (L) = min{e, (lOeMiL^)" 1 }. (5.17) 

For z G B> € (l)(zo) and n G Q', we then have a n (z)^ < r/4 and hence Z' n (TL, z) = Z ti (Tl, z), 
implying in particular that 

~ Q (z) = Z h i g (z)+ Z mO?L,z). (5.18) 

me.S\S' 

Note that this implies, in particular, that Sq(-) is analytic in D e (L) (zq). 

Our next goal is to prove a suitable bound on the right hand side of (5.18). By Lemma 4. 10, the 
first term contributes no more than 2L d C{z) Ld e- TL / i , provided r > 4c +16andL is so large that 
5L d e -rL/4 < log 2 on the other hand, since z G B e(L ) {z ) implies that that a m (z) > r/(16L) 
for all m Q! , the bound (4.74) implies that each Z m (Ti, z) on the right hand side of (5.18) 
contributes less than 2((z) L<I e~ Tljd X / 32 once L is so large that AL d e~ TL l A < log 2. By putting 
all of these bounds together and using that Q{z) Ld < ((z ) Ld e Ml \ z ~ z °\ Ld < e l /^C,(z Q ) Ld by 
the bound (4.63) and our definition of e^ L \ we get that 

\S Q (z)\ < 5\S\L d C(z ) Ld e~ rLd - 1 ^ 2 (5.19) 

whenever z G D e ( L )(zo) and L is so large that L > r/2 and 5L d e~ rL / 4 < log 2. Increasing L if 
necessary to guarantee that e( L ) = (10eMiL rf ) _1 and applying Cauchy's theorem to bound the 
derivatives of Eq(z), we thus get 

dl~ Q (z) < £!(10 e M 1 ) £ 5|5|L^+ 1 )c(^o) Ld e" rL ^ 1/32 (5.20) 

z=zo 

provided L > Lq, where Lq = Lo(d, M, r, e) is chosen in such a way that for L > Lq, we have 
L > r/2, 5L d e- rL / 4 < log 2 and (WeMiL 11 )' 1 < e. Since z G W k /l{Q) was arbitrary and 
|<S| = YtmeK 1m, this proves the desired bound (2.20) with C = lOeMi = 10e(l + AM). □ 
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